Time-dependent angularly averaged inverse transport 

extended version 



o 
o 

(N 
>v 



Guillaume Bal and Alexandre Jollivet 
May 7, 2009 

Abstract 



This paper concerns the reconstruction of the absorption and scattering parameters in 
a time-dependent linear transport equation from knowledge of angularly averaged mea- 
surements performed at the boundary of a domain of interest. We show that the absorption 
coefficient and the spatial component of the scattering coefficient are uniquely determined 
i-C . by such measurements. We obtain stability results on the reconstruction of the absorption 

and scattering parameters with respect to the measured albedo operator. The stability 
results are obtained by a precise decomposition of the measurements into components 
with different singular behavior in the time domain. 



(M '■ 1 Introduction 

> : 

■ Inverse transport theory has many applications in e.g. medical and geophysical imaging. It 
"st" ■ consists of reconstructing optical parameters in a domain of interest from measurements of 
^} . the transport solution at the boundary of that domain. The optical parameters are the total 
^ ! absorption (extinction) parameter o{x) and the scattering parameter k(x, v', v), which measures 
Q\ \ the probability of a particle at position i 6 I C I" to scatter from direction v' G S n_1 to 
O ■ direction v G § n_1 , where § n_1 is the unit sphere in M n . 

The domain of interest is probed as follows. A known flux of particles enters the domain 
^ . and the flux of outgoing particles is measured at the domain's boundary. Several inverse theories 
may then be envisioned based on available data. In this paper, we assume availability of time 
dependent measurements that are angularly averaged. Also the source term used to probe 
the domain is not resolved angularly in order to e.g. save time in the acquisition of data. 
More precisely, the incoming density of particles 4>(t, x, v) as a function of time t, at position 
x G dX at the boundary of the domain of interest, and for incoming directions v, is of the 
form (fis(t,x,v) = <p(t,x)S(x,v), where <f)(t,x) is arbitrary but S(x,v) is fixed. This paper is 
concerned with the reconstruction of the optical parameters from such measurements. We show 
that the attenuation coefficient is uniquely determined and that the spatial structure of the 
scattering coefficient can be reconstructed provided that scattering vanishes in the vicinity of 
the domain's boundary (except in dimension n = 2 and when A is a disc, where our theory 
does not require k to vanish in the vicinity of OX). For instance, when k(x, v', v) = ko(x)g(v', v) 
with g(v',v) known a priori, then ko(x) is uniquely determined by the measurements. Similar 
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results were announced in pQ when measurements are available in the modulation frequency 
variable, which is the dual (Fourier) variable to the time variable. 

Several other regimes have been considered in the literature. The uniqueness of the recon- 
struction of the optical parameters from knowledge of angularly resolved measurements both 
in the time-dependent and time-independent settings was proved in [HI 10] ; see also [T7j for a 
review. Stability in the time-independent case has been analyzed in dimension n = 2, 3 under 
smallness assumptions for the optical parameters in O [T5] and in dimension n = 2 in [18J. 
Stability results in the presence of full, angularly resolved, measurements have been obtained 
in [HUH [19]. The intermediate case of angularly averaged measurements with angularly resolved 
sources was considered in [13j. The lack of stability of the reconstruction in the time indepen- 
dent setting with angularly averaged measurements and isotropic sources is treated in [6]. See 
also [2] for a recent review of results in inverse transport theory. 

The rest of the paper is structured as follows. Section [5] recalls known results on the 
transport equation and the decomposition of the albedo operator. In section [3] we define and 
decompose the averaged albedo operator (Proposition 13.11) and we study its distributional ker- 
nel (Theorems I3.2fl3~5l . Our main results on uniqueness and stability are presented in section 
141 (Theorems I4TTH4.21 Theorems 14. #14751 and Corollary 14.61) . We show that the absorption co- 
efficient and the spatial structure of the scattering coefficient (the phase function describing 
scattering from v to v' has to be known in advance) can be reconstructed stably from angularly 
averaged time dependent data. The reconstruction of the scattering coefficient requires invertion 
of a weighted Radon transform in the general case. In the specific case of a spherical geometry 
(measurements are performed at the boundary of a sphere), then the scattering coefficient may 
be obtained by inverting a classical Radon transform. In section [5] we prove Theorems I3.4H3.51 
In section [H] we prove Theorems 14. H - T4~2"l Theorem 14.51 and Theorem 14.41 f)4.13p . In section [7] we 
prove Theorem 13 .21 In section [8] we prove Theorem 13.31 In section [9] we prove Lemmas I8.1H8.4I 
that are used in section [HJ In section [10] we prove Proposition 13.11 

The derivation of the results is fairly technical and is based on a careful analysis of the 
temporal behavior of the decomposition of the albedo operator into components that are multi- 
linear in the scattering coefficient. Our results are based on showing that the ballistic and single 
scattering components may be separated from the rest of the data. These two components are 
then used to obtain our uniqueness and stability results. It turns out that the structure of single 
scattering is different depending on whether k vanishes on dX or not. When k does not vanish 
on dX, the main singularities of the single scattering component do not allow us to "see inside" 
the domain as they only depend on values of k at the domain's boundary in dimension n > 3. 
The singular structure of single scattering and the resulting stability estimates are presented 
in detail when both k vanishes and does not vanish on dX. 

This is the extended version of a submitted paper [5]. 

2 The forward problem and albedo operator 
2.1 The linear Boltzmann transport equation 

We now introduce notation and recall some known results on the linear transport equation. 
Let X be a bounded open subset of M. n , n > 2, with a C 1 boundary dX. Let z/(x) denote the 
outward normal unit vector to dX at x G dX. Let T± = {(x, v) G dX x § n_1 | ± v(x) ■ v > 0} 
be the sets of incoming and outgoing conditions. For (x, v) G X x S n_1 we define t±(x, v) and 
t(x, v) by t±(x,v) := inf{s G (0, +oo) | x ± sv ^ X} and t(x,v) := r_(x,t>) +t + (x,v). For 
x G dX we define := {v G S n ~ x | ± v(x) ■ v > 0}. 
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Consider a : X x §™ 1 — > R and : X x §™ 1 x §™ 1 — > M. two nonnegative measurable 
functions. We assume that (a, k) is admissible when 



< a G L°°(X x S™" 1 ), 

< t/, .) G L^S* -1 ) for a.e. (x, v') G X x S n_1 , 
a p (x, u') = / k(x, v', v)dv belongs to L°°(X x S"" 1 ). 

Let T > i] > 0. We consider the following linear Boltzmann transport equation 
du 

— (t, x, v) + v ■ V x u(t, x, v) + a(x, v)u(t, x, v) 

= / k{x,v',v)u{t,x,v')dv', (t,x,v) G (0,T) x X x S n_1 , 

M|( ,T)xr_ (*, a;, = (pit, x, v), 



(2.1) 



(2.2) 



u(0,x,t>) = 0, G X x § 



n— 1 



where G L 1 ((0 , T) , L 1 (T _ , d£)) and supp</> C [0,7/]. Here, d£(x, v) = \v ■ u(x)\dvdfj,(x), where 
dfi is the surface measure on dX and dv is the surface measure on S n_1 . In other words, we 
assume that the initial condition is concentrated in the ^-vicinity of t = and measurements 
are performed for time T, which we will choose sufficiently large so that particles have the time 
to travel through X and be measured. 

2.2 Semigroups and unbounded operators 

We introduce the following space 

Z := {f eL\X xS n -') rVJe^flxS"- 1 )}, (2.3) 

Lipfxs™- 1 ) + \\v ■ VxfW^iXxS™- 1 )', (2.4) 



where v ■ V x is understood in the distributional sense. 

It is known (see [S]) that the trace map 7_ from C 1 (X x S™" 1 ) to C(r_) defined by 

7-(/) = /|r_ (2.5) 

extends to a continuous operator from Z onto L 1 (T_, r + (x, v)d£(x, v)) and admits a continuous 
lifting. Note that L 1 (r_,c/^) is a subset of the spaces L x (r_, t + (x, v )d£(x, v )). 
We introduce the following notation 

AJ = -af, A 2 f = [ k(x, v', v)f(x, v')dv'. (2.6) 



As (a, k) is admissible, the operators A\ and A 2 are bounded operators in L 1 (X x§" 1 ). 
Consider the following unbounded operators 

71/ = -v ■ VJ + A t f, D{T X ) = {feZ\ /, r _ = 0}, (2.7) 
Tf = T x f + A 2 f, D(T) = D^). (2.8) 

It is known that the unbounded operators T\ and T are generators of strongly continuous 
semigroups in L l (X x S n_1 ) U\(t), U(t) respectively (see for example [HI Proposition 2 pp 
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226]). In addition U\(t) and U(t) preserve the cone of positive functions, and U\{t) is given 
explicitly by the following formula 

U x {t)f = e'^^- s ^) ds f( x -tv,v)e(x,x-tv), for a.e. (x,v) e X x S™" 1 , (2.9) 
for /el^IxS"- 1 ), where 

,A _ / 1 if x + P(2/ ~ x ) e X for a11 P e (°> !]> / 91 m 

UlX ' yJ ~\ otherwise, 1 UJ 

for (x,y) ei"xK" 

We recall the Dyson-Phillips formula 

+0O 



U(t)=J2 H ™(t) (2.H) 



m=0 

for t > 0, where 



F m (t) := J Ui(t- st- s m )A 2 Ui(s l ) . . . A 2 U 1 (s m )ds 1 . . .ds m , m > 1, (2.12) 



s 1 >0,...,s m >0 
sx + --- + s m<* 

ft 



# ro (t) = f H m -i(t-s)A 2 U 1 (s)ds, m>l, (2.13) 
#o(*) := ^(t). (2.14) 



2.3 Trace results 

We introduce the following space 



W := {u E L l {{0,T) x X x S"- 1 ) | ^ + o.V I jHGL 1 ((0,T)xIxr 1 )}, (2.15) 



|m||w := ll^llii^o.TJxXxS"- 1 ) + 



9 ^ ^ 



LUtO.TJxXxS"- 1 ) 



(2.16) 



where and u • V x are understood in the distributional sense. 

It is known (see [3 IE]) that the trace map 7_ (respectively 7+) from C 1 ([0,T] x X x § n_1 ) 
to C(X x S"- 1 ) x C((0, T) x r±) defined by 

7-0) = 0(0, .),^|( ,T)xr_) (respectively 7+^) = (ip(T, .), V>|(o,T)xr+) (2.17) 

extends to a continuous operator from W onto L X (X x § n_1 , r + (x, v)dxdv ) x L x ((0, T) x 
r_,min(T - t,r + (x,v))dtd£(x,v)) (respectively L X {X x S n_1 , r_(ar, v)dxdv) x ^(((^T) x T+, 
min(t, t_(x, v))dtd£(x, v ))). In addition 7^. admits a continuous lifting. Note that L X (X x § n_1 ) 
is a subset of L l (X x § n -\ T+ (x, v)dxdv). Note also that L 1 ((0, T) x T_,dtd£) (resp. L 1 ((0, T) x 
r + , dtd£)) is a subset of -^ 1 ((0, T) xT_, min(T— t, r + (x, v))dtd£(x, v)) (resp. ^ 1 ((0, T) xT + , min(t, 
r_(x, v))dtd^{x, v ))). 

We now introduce the space 

VV := {u e W I 7_(«) G L X (X x S"" 1 ) x L\(0,T) x r_,dtd£)}. (2-18) 

We recall the following trace results (owed to jH E] in a more general setting). 
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Lemma 2.1. The following equality is valid 

W = { u e W | 7+(u) 6 x S"" 1 ) x L x ((0, T) x r+, dtd£)}- ( 2 - 19 ) 
In addition the trace maps 

7 ± : VV — > L 1 (Xx§ n_1 ) xL 1 ((0, T) xT-t, are continuous, onto, and admit continuous lifting. 

(2.20) 

2.4 Solution to equation ( 12.21 ) 

We identify the space I/ 1 ((0, r), L x (r±, d£)) with the space L 1 ((0,r) x r.|-,dtd£) for any r > 0. 
We extend by on M. outside the interval (0,rj) any function G -^ 1 ((0, rj), L 1 (T_, d£)). 

Let G L 1 ((0,r/),L 1 (r_,c/O)- Th en we consider the lifting G-(t)<j> G VV of (0,0) defined 

by 

G-(t)<j>(x,v) := e-f°~ lX ' V)eT{x - av ' v)ds (l)-(t-T_(x,v),x-T_(x,v)v,v), for a.e. (t,x,v) G (0, T) xXxS n_1 . 

(2.21) 

Note that G_ (.)0 is a solution in the distributional sense of the equation ( + v ■ V x )u + cru — 
in (0, T) x X x S™- 1 and 

||G_(.)0||w < (1 + ||cr||oo)||G'_(.)0|| L i (( o i T)xXx§™-i) < (1 + IM|«O T ll0-IU 1 ((o,. J )xr_,<te(e)- (2.22) 

To prove this two latter statements, one can use the following change of variables (see [ID])- 

Lemma 2.2. We have 

r r rr±(x,v) 

/ f(x,v)dxdv= / / f(x±tv)dtd^(x,v), (2.23) 

JxxS"- 1 Jt t Jo 

for f G L X (X x V). 

From QZZZD we obtain that the map i : ^((0, rj), L 1 ^-, df)) -> VV defined by 

z(0) = G_(.)0, 0GL 1 ((O,r ? ),L 1 (r_,dO), (2-24) 

is continuous. 

The following result holds (see [TT[ Theorem 3 p. 229]). 
Lemma 2.3. JTie equation (12. 2p admits a unique solution u in W which is given by 

u(t) = G4t)(f)+ [ U(t- s)A 2 G_(s)(j)ds. (2.25) 

where U(t) is the strongly continuous semigroup in L 1 (X x §>™ _1 ) introduced in subsection \2.2[ 

Using (12.251) and the Dyson-Phillips expansion (12.111) we obtain that the solution u of (12.21) 

may be decomposed as 

u(t) = G-(t)<p +J2 Hm ^ ~ s ) A 2G-(s)(j)ds, (2.26) 

m =0 J-°° 

for t > and G L 1 ((0,?7) x dX , dtdfi(x)) . The first term in the above series G_(t)0 is the 
ballistic part of u(t) while the term corresponding to m > 1 is m-linear in the scattering kernel 
k. The term corresponding to m = 1 is the single scattering term. 

From (12.241) . Lemma [2.31 and (12.201) . we also obtain the existence of the albedo operator. 
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Lemma 2.4. The albedo operator A given by the formula 

A(f) = M|(o,T)xr+, for 4> £ -^ 1 ((0, rj), L 1 (T_, d£)) where u is given by (12.251) . (2.27) 
is well-defined and is a bounded operator from -^ 1 ((0, 77), L 1 (T_, d£)) to -^ 1 ((0, T), L 1 (T + , d£)). 

We refer the reader to [9] for the reconstruction of the optical parameters when the full 
albedo operator is known. We assume here that only partial knowledge of the albedo operator 
is available from measurements. 

3 The operator Agyy and its distributional kernel 
3.1 Angularly averaged measurements 

We now define more precisely the type of measurements we consider in this paper. The direc- 
tional behavior of the source term is determined by a fixed function S(x, v), which is bounded 
and continuous on T_. We assume that the incoming conditions have the following structure 

4> s (t',x',v') = S(x',v')4>(t',x'), t' e (0,77), (x',v') e r_, (3.1) 

where <f>(t, x) is an arbitrary function in -^ 1 ((0, 77) x dX). We model the detectors by the kernel 
W(x,v), which we assume is a continuous and bounded function on T + . The available mea- 
surements are therefore modeled by the availability of the averaged albedo operator As,w from 
^((0,77) x dX,dtd{i(x)) to ^((O,! 7 ) x dX,dtdfi(x)) and defined by 

A StW <f>{t,x) = / A((f)s)(t,x,v)W(x,v)(u(x) -v)dv, for a.e. (t,x) G (0,T) x dX. (3.2) 



n-l 



The functions S and W are fixed throughout the paper. The case W = 1 corresponds to 
measurements of the current of exiting particles at the domain's boundary. 

The decomposition of the transport solution (I2.26P translates into a similar decomposition 
of the albedo operator of the form 



+00 



A SiW <j>(t, x) = ^2 A ™^w<P(t, x), (3.3) 

m=0 

for (t, x) G (0, T) x dX, where we have defined 

A 0>s ,w<l>(t, x) = / (u(x) ■ v)W(x,v) (G ? -(-)0s)|( O) T)xr + {t,x,v)dv, (3.4) 



n-l 



A m> s,w<f>(t,x) = {u{x) ■ v) W(x, v) I / H m -i(t - s)A 2 G-(s)(f)sds ) (t,x,v)dv, (3.5) 

V./-00 / |(o,T)xr+ 

for a.e. (t,x) G (0, T) x dX where (fi s is defined by (13.11) . The kernels of the operators A m)S> w 
can be written explicitly. 
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3.2 Distributional kernel of the operators A mj s,w 
Consider the nonnegative measurable E from dX x dX — > K defined by 

e -Jo i=i-= 2 i' i*i -* 2 i s if xi + p(x 2 - x x ) G X for all p G (0, 1), 

otherwise, 

for a.e. (xx, x 2 ) G dX x dX. For m > 3, we also define the nonnegative measurable real function 
E(xi, . . . , x m ) by the formula 

=e ^ jo a{x * 'i-i-mlM^-i+ii^e^^^O^eCxi.Xi+x), (3.7) 

for a.e. (xi,...,x m ) e Six (M™)" 1 - 2 x <9X, where 6 is defined by (l21Dj) . The function 
E(xi, . . . ,x m ) measures the total attenuation along the broken path (xi, . . . ,x m ) G dX x 
Rm -2 x dx provic ied (x 2 , • • • , x m _i) G X m - 2 . 

For m G N, m > 1 and for any subset U of W 1 we denote by xt/ the characteristic function 
from R m to K. defined by Xu{y) — 1 when y G U and Xu(y) — otherwise. Using (I3.4p - fl3.5l) . 
f !2.2ip and fl2.13l) - fl2.14l) we then obtain the following result on the structure of the kernels of 
the albedo operator. 

Proposition 3.1. We have 

A m ,s,w(<P)(t,x) = / J m (t-t' \x,x')<j>(t',x')dt'dn(x'), (3.8) 

J(0,ri)xdX 

form > and for a.e. (t,x) G (0, T) x <9X ; where 
l0 (r,x,x') := ,^^2 1 -tOK^O -HU^^t- \x-x'\), (3.9) 

7i(t, x, x') := X(o +oo)( T " — W — x\) / (y(x) ■ v )W(x, v ) [E(x, x — sv , x')k(x — sv, v ', v ) 

2 n - 2 (Y — (x — x') • f ) n ~ 3 

X X(Q,T-(x,v))(s)S(xr,v')\u(a?)-V , \]i_ x _ x , T a_ k _«>|2 |^_^_^,|2n-4 du » ( 3 - 10 ) 

for (r, x, x') G Mx dX x <9X and where 7 m /or m >2 admits a similar, more complex, expression 
given in Section^ (see fl8~T2l) - fl8~T3l ). 

Because the above formulas are central in our uniqueness and stability results, we briefly 
present their derivation and refer the reader to Section [10] for the rest of the proof of Proposition 
O 

Derivation of (13.91) and (13.101) . From (13.41) and the definition of G_, we obtain Aq^^v^^iX) = 
L n -i(v(x)-v)W(x, v)E(x, x — r_(x, v)v)S(x— r_(x, v)v, v)(p(t—r_(x, v), x— r_(x, v)v)dv, {t, x) G 

(0,T) x dX and for G L 1 ((0,?7) x dX). Therefore, performing the change of variables "x"' = 
x — t(x, f )u {dv = ^}^n-i dn(x') and r(x, v) = \x — x'\), we obtain (13.91) . 

From the definition of A 2 and G_ we note that A 2 G^(s)(ps(z, w) '■= Ln-i k(z, v', w)E(z, z — 
t_(z, v')v')S{z — T-(z, v')v', v')(f>(s — T-(z, v'), z — r_(z, v' ) v' ) dv' , for a.e. (z, w) G X x § n_1 and 
for G L 1 ((0,n) x <9X). Performing the change of variables "x' = z — r_(2:, ?/)?/" , we obtain 
the equality (A 2 G_(s)0 s ) (z,to) = f gx [k(z,v' ,w)S{x' ,v')\v{x') ■ v'\] v>= z - x > r§^k <fr(s - \z - 

\z — x' I ' ' 

x'|, x')dfi(x'), for a.e. (z, w) G X x S n_1 and 4> G L 1 ((0,n) x 9X). Using also the definition of 
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Ai t s,w ( see (13.51) for m = 1) we obtain the following equality for any G ^((O,^) x dX) and 
for a.e. if, x) G (0, T) x <9X 

A hS ,w{<P){t,x) = / / [fc(g- (t - sKt/^fffx'^OKx') 'i/IL rMa! ( v ( x ) ' v ) 



x,+ 



V j j- r j- 

" -it— S)V — X 



x - — — ^ — — ^ — \®{x, x — (f — s)v)(j)(s — \x — if — s)u — x'l, v)dfi(x')dsdv. (3.11) 

|x — [t — s)v — x'\ n ~ L 



Then performing the changes of variables "s" = t — s and "i"' = i — s — \x — sf — x'| (s = 

(t-t') 2 -\x~x'\ 2 dt' 2((t-t')-(x-x')-v) 2 \ i , • /nTTTTh n 

2{t-J-v{x-J)V ITS = \x-J-lt-t')v\ 2 )» We ° btam CT)' D 

To simplify notation, we define the multiple scattering kernels 

oo 

r fc = ^ 7m - (3.12) 



m=k 



3.3 Regularity of the albedo kernels 

The reconstruction of the optical parameters is based on an analysis of the behavior in time of 
the kernels of the albedo operator. Our first result in this direction is the following. 

Theorem 3.2. Assume that k G L°°(X x S™" 1 x S^ 1 ). Then the following holds: 

^t 2 - \x-x'\ 2 k{t,x,x') G L°°((0,T) x ^ x dX) when n = 2; (3.13) 
4y7i {t, x, x') G L°°{{0, T) x 3X x dX) when n = 3; (3.14) 



T\X — X 



\ n I r+\x-x'\ 
V r— |ai— x'\ 

t\x- x'l^'^i^x,^) E L°°((0,T) x OX x OX) when n > A. (3.15) 

In addition, assume that k G L°°(X x § ra_1 x S n_1 ) and that there exists 5 > st/c/i £/ia£ 
supp/c C {y g X | inf^x |# — 2/| > «-e., £/ie scattering coefficient vanishes in the vicinity 
of dX . Then, the following holds 

(r- \x-x'\) — 7l (r,x,x / ) G L°°((0,T) x dX x dX) when n > 2. (3.16) 



Theorem l3.2l is proved in Section[71 The results (13. 14ft and (13.15p of Theorem [32] correspond 
to singularities of the single scattering contribution that depend on the values of k on dX. The 
above theorem shows that the structure of the single scattering coefficient is quite different 
depending on whether k vanishes on dX or not. 

The following result describes some regularity properties of the multiple scattering. It is 
because multiple scattering is more regular than single scattering, in an appropriate sense, that 
we can reconstruct the scattering coefficient in a stable manner. 

Theorem 3.3. Assume that k G L°°(X x S^ 1 x S"" 1 ). Then the following holds: 

T 2 (r, x, x) G L°°((0, T)xdXx dX), when n = 2; (3.17) 
T \ X - X '\ T ^ T ^ X ') _ e L oo ((Q; y) x qx x dX), when n = 3; (3.18) 



1 + ln 

Tlx — X 



T+\x— X'\ 



T — \X — X 



T—\x—x'\ 

■ n-2 

—T 2 (t,x,x') G L°°((0,T) x dX x dX), when n > A. (3.19) 



S 



In addition, assume that 6 L°°(I x S"- 1 x S™- 1 ) and that there exists 5 > such that 
supp/c C {y e X | inf^gdx \x — y\ > 5}. Then the following holds: 

(r-lx-^l)' 1 ^l + ln ^ + j^~^j ^ T 2 (t,x,x') G L°°((0,T) x OX x dX), when n = 3; 

(3.20) 

1 — 71 

(T-\x-x'\) — T 2 (T,x,x , )eL oo ({0,T)xdXxdX), whenn>A. (3.21) 

Theorem 13.31 is proved in Section [HJ These results quantify how "smoother" multiple scat- 
tering is compared to the single scattering contribution considered in Theorem 13.21 

3.4 Asymptotics of the single scattering term 

In this subsection we assume that X is also convex. We give limits for the single scattering 
term in two configurations given by: 



.1 J\ A. L 

the nonnegative function k is continuous on X x § n_1 x 



the nonnegative function a is bounded and continuous on X x §" 

1 gn- 1. 



(3.22) 



or 



there exists a convex open subset Y C X with C 1 boundary such that a(x,v) = 

for (x,v) G (X\Y) x S n_1 and the nonnegative function a is bounded and 

continuous on Y x §>™ _1 ; and there exists a convex open subset 

with C 1 boundary such that 5 := in£( x ,z)edXxZ \x — z\ > and k(x, v',v) = 

for (x,v',v) G (X\Z) x § n_1 x 8 n_1 , and the nonnegative function k 

is bounded and continuous on Z x S n_1 x S™ -1 , 



(3.23) 



When either (13.221) or (13.231) is satisfied, we want to analyze the behavior of the function 
7i(t, x, x') given by the right hand side of (13.101) for all (r, x, x') G K x <9X x <9X. We need to 
introduce some notation. Let i?o : S n_1 x X — > K be the function defined by 

0„(v,a;) = (r_(x,w)r + (x,t;))-^, (v,x) G § n_1 x X, (3.24) 

and consider the weighted X-ray transform P^ defined by 

PT+ (x,v) 

Pt f(v,x)= tf {v,tv + x)f(tv + x)dt, (3.25) 

J T_ (x,v) 

for a.e. (v,x) G S n_1 x dX and / G £ 2 (X, sup,, g§ n-i f?o(i;, x)dx). The first result analyzes the 
behavior of 71 under hypothesis (13.221) . 

Theorem 3.4. Assume that the open subset X ofM. n with C l boundary is convex. Let (x, x' ) G 
dX 2 be such that x + s(x — x' ) G X for some s G (0,1). Set t> = Z x °\ an d to = \x — x' \. Then 
under condition (13.221) . we have the following results. When n = 2, then 

, 1 \Z2W(x,v )S(x , ,v )(v(x) ■ v )\is(x ) ■ v \E(x,x' ) 

7il r ) x i x o) - 



V r ~to (3.26) 
xP# k V0 (v ,x) +0 (—===), asr^t^ 

\\/T — in/ 
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where P$ is defined by (13.251) and k vo (y) := k(y,v ,v ) fory G X. 
When n = 3, then 



7i(r,x,x' ) 
When n > 4, then 



ln( 



1 



7T . 



/7I W(x,vo)S(x' ,vo)(v(x) ■v )\p{x' Q ) ■ v \E(x,x' ) 

T-t t z Q 



x (k(x, v , v Q ) + k(x' ,v ,v )) +o( ln( — —)), as r -»• 

V r — t / 



t E(X) Xq) 



S(x' ,v )\p(x' o ) ■ Vq 



W(x,v)(v(x) ■ v)k(x, Vq, V ) 



(3.27) 



n-l 



+W(x,V )(l/(x) ■ Vq) 

+o(l), ^ . 



n— 1 
/ _ 
0' 



1 - W 
1 — W ' ■ f o 



dv 



dv' 



(3.28) 



Theorem 13.41 is proved in Section [51 Note that 71 depends on the value of k on dX in 
dimension n > 3. Under hypothesis (13.231) . i.e., when the scattering coefficient vanishes in the 
vicinity of where measurements are collected, we have the quite different behavior: 

Theorem 3.5. Assume that the open subset X of R n with C 1 boundary is also convex and 
assume that condition (I3.23P is fulfilled. Let (x,x' ) G dX 2 be such that x' Q + s(x — x' ) G Z for 
some s G (0, 1). Set v = ZZ^n an d ^0 = \ x ~ x 'o\- Then we have the following. 

When n = 2, then (EOSj) still holds. 

When n > 3, then 

-yx(r,x,x' ) = (r-t )^(2t )^Vo\ n ^ n - 2 )S(x' ,v )W(x,v )\u(x^ 

n — 3 

xE(x,x' o )P# o k Vo (v 0} x) + o((r - t )~s~), as r -> tj, 
where P$ is defined by (13.251) and k VQ (y) := k(y, vq, vq) for y G X . 



(3.29) 



Theorem 13.51 is proved in Section Theorem 13.51 may remain valid under different condi- 
tions from those stated in (13. 23ft . For instance, when a is bounded and continuous on X and 
k is continuous on X x S n_1 x S n_1 and k(x, ., .) decays sufficiently rapidly as x get closer and 
closer to the boundary dX for any x G X, then the same asymptotics of 71 holds. 



4 Uniqueness and stability results 

We denote by 7 := T = X]m=oT m ^he distributional kernel of A S:W . Then 7 — 70 = I\ 
denotes the distributional kernel of the multiple scattering of A$ } w- F° r the res t of the paper, 
we assume that the duration of measurement T > diam(X) := sup( 2 . y - )gX 2 \x — y\ so that 
the singularities of the ballistic and single scattering contributions are indeed captured by the 
available measurements. 

Let (a, k) be a pair of absorption and scattering coefficients that also satisfy (12. ip . We 
denote by a superscript" any object (such as the albedo operator A or the distributional kernels 
7 and To) associated to (a, k). Moreover if (ex, k) satisfies (13.231) for some (Y, Z) and (a, k) also 
satisfies (13.231) for some (Y,Z), then we always make the additional assumption Y = Y and 

Z = Z. Let ll-II^T := \\-\\c(L^((0,v)xdX)),L^((0,T)xdX))- 
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4.1 Stability estimates under condition ( 13.221 ) or ( 13.231 ) 

Theorem 4.1. Assume that the open subset X of M n with C 1 boundary is also convex. Let 
(a,k) and (a,k) satisfy condition ( 13.22)) . Let x' G dX. Then we have: 



dX 



< \\A s ,w ~ A s ,w\\r,,T- 

t =\x-sc'\ 



(4 ' 1} 

Let x G <9X 6e suc/i that px' + (l—p)x G X /or some p G (0, 1). Set i>o = ]fzf?T anc ^ = | x — o: 1 . 
When n = 2, we have 



W(x,v )S(x' ,v )(u(x) ■v )\u(x' ) -v \ E(x,x' )P A) k Vo (v ,x) - E(x, x' )P# k Vo (v , x] 
1 y / T ^-\z-z'\^r 1 -r 1 )(T,z,z') , 



< - 

- 2 



(4.2) 



where \\ ■ \\ L oo := \\ ■ \\ L ™((o,T)xdXxdX), P$ is defined by ()3.25p and k Vo (y) := k(y,v ,v ) for 
y G X (k vo is defined similarly). 
When n = 3, then 



E(x, x' )(k(x, v , v ) + k(x' Q , v , v )) - E(x, x' Q ) (k(x, v , v ) + k(x' , v , v )) 



xW(x,v )S(x' ,v )(v(x) -v )\v(x! ) -v \ < - 

7T 



T\Z — Z 



r— \z — z'\ 



i^-roew') 



.(4.3) 



When n > 4, then 
S(a? ,v )\u(aQ ■ v \ 

+W(x,v )(v(x) -v ) 



W(x, v)(u(x) ■ v) 

-i 1 — v ■ V o 

f 

S(x' ,v')\u(x' )-v' 

n-1 1 — V' ■ VQ 



E(x, x' )k(x, v o, v) — E(x, x' )k(x, vq, v) ) dv 



< 



T\z-z'\ n - 2 (ri-ri)(T,z, z') 



E(x, x' )k(x' , v v ) — E{x, x' Q )k(x' , v', v ) J dv' 



(4.4) 



Theorem 14.11 is proved in Section [71 It shows that the spatial structure of k may be stably 
reconstructed at the domain's boundary. More interesting is the following theorem, which pro- 
vides some stability of the reconstruction of the scattering coefficient when it vanishes in the 
vicinity of the boundary dX. 

Theorem 4.2. Assume that the open subset X o/R n with C 1 boundary is also convex. Assume 
also that inf(a;/y) e r_ S(x', v') > and inf( X) „) e r + W(x, v) > 0. Let (a,k) and (a,k) satisfy 
condition (13.231) . Let x' G dX. Then there exist constants G\ = Ci(S,W,X,Y) and C 2 = 
C 2 (S } W } X } Z) such that 



n-l 

I _ 
0' 



\E - E\(x' + t + (x' ,vo)vo,x' )\v(x' ) ■ v \dv < Ci\\A s , w - A Si w\\ v ,t, (4.5) 



E(x, x' )P# o k v/o (v' , x' ) - E(x, x' )P# o k v > o (v , x' ) 



[r-\z-z'\)-* t (r 1 -r 1 )(T,z,z') 



L°° 

(4.6) 
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for x G dX such that px' + (1 — p)x G Z for some p G (0, 1) where v' = jfzfrj; P&a ^ s defined 
in (13.251) . and k v > o (y) := k(y,v' Q ,v' Q ) fory G X (k v > is defined similarly). 

Theorem 14.21 which is one of the main results of this paper, is proved in Section [HJ 
4.2 The case when X is a ball of R n 

When X is an open Euclidean ball of M n , which is important from the practical point of 
view in medical imaging as it is relatively straightforward to place sources and detectors 
on a sphere, we are able to invert the weighted X-ray transform P# f, f{x,v) := f(x) G 
L 2 (X, sup^ggn-i $o(v, x)dx) using the classical inverse X-ray transform (inverse Radon trans- 
form in dimension n — 2). In the next subsection we shall consider a larger class of domains 
X, which requires one to solve more complex weighted X-ray transforms. 

Up to rescaling, we assume X = B n (0, 1), the ball in W 1 centered at of radius 1. Consider 
the X-ray transform P defined by 

PT + (X,V) 

Pf(v, x) = / f(sv + x)ds for a.e. (v, x) G S™" 1 x dX, (4.7) 

J T- (x,v) 

for / G L 2 (X) (we extend / by outside X). We have the following Proposition 14.31 
Proposition 4.3. When X = B n (0, 1) we have 

PtJ(v,x)=P( 6 f)(v,x), for a.e. (v, x) G S^ 1 x dX, (4.8) 

for f G L 2 (X,swp v(ES n-i'do{v,x)dx) where g(y) := (1 - |y| 2 )"V ; y E X. 

Proof of Proposition \4-3[ It is easy to see that 

T±{tv + qv ± ,v) = y/l-q 2 =ft, (4.9) 
■&o(v,x) = (l-q 2 -t 2 )-^ = (1- \x\ 2 )-^, (4.10) 

for {t, q) G M 2 , t 2 + q 2 < 1 and for (v , v 1 ) G § n_1 x S n-1 , v ■ v 1 = 0, where x — tv + qv 1 (we 
remind that i?o is defined by (13. 24ft ). Then Prop osit ion 14 . 3 1 follows from the definition ( 13.25D . □ 

Assume that (cr, k) satisfies condition (I3.22p when n = 2 or (I3.23P when n > 2. Assume 
also that k(x,v,v') = k (x)g(v,v') for a.e. (x,v,v f ) G X x S n_1 x S n_1 where g is a given 
continuous function on S n_1 x S™^ 1 , inf„ 6 §n-i g(v,v) > 0, and where k G L°°(X). Then from 
the decomposition of the angularly averaged albedo operator As t w (Proposition 13. ip and from 
Theorems 13.21 13.31 13.41 and 13.51 and from Proposition 14.31 and methods of reconstruction of 
a function from its X-ray transform, it follows that (cr, ho) can be reconstructed from the 
asymptotic expansion in time of A s> w provided that a = a(x) and inf(yy) gr S(x',v') > 
and inf( :E) „) e r + W(x, v) > 0. In addition we have the following stability estimates. 

Theorem 4.4. Assume X = B n (0, 1) and inf(yy) 6r _ S(x',v') > and inf( Xjl ,) gr+ W(x, v) > 0. 
Let (o~,k) and (a,k) satisfy either condition (I3.23P or (13.221) . Assume that a, a do not depend 
on the velocity variable (cr(x, v) = a(x) ) and suppaUsuppa C Y , where Y C X is a convex open 
subset ofW 1 with C 1 boundary, and let M = max(||c||£oorn, 1 1 5" 1 1 ) . Assume k(x,v,v') = 
k (x)g(v,v') and k(x,v,v') = k (x)g(v , v '), g(v,v) > 0, for (x,v,v f ) G X x S n ~ 1 x § n_1 where 
g is an a priori known continuous function on § rt_1 x § ra_1 . 
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Then there exists C3 = Cs(S, W, X, Y, M) such that 



a — a 



1 

~ < C3\W — <5'lll°°(Y)||As , ,W' r — As^wHI^T- 



(4.11) 



When n > 2 and (a, k) satisfies (I3.23p . there exists Ca,\ = C^i(S, W, X, Y, Z, M, g) such that 

(4.12) 

1 

2 



g(k — k )\\ H _i^ < C4 t i\\k — k \\So y\\k \\oo\\As,w ~ ^wW^t 



(T-lz-^T^i-fiXT,*,^ 



When n = 2 and (a, A;) satisfies (13.221) . there exists C^ 2 = C^ 2 (S, W, X, M, g) such that 

I ^0 II 00 1| A s, w ~ A s, w II ti,t (4-13) 



\g(k - k )\\ H _i < C 4)2 \\k - k 



3 

1 4 

OO 



Iz-zf^-roo-,*,*') 



Theorem 14.41 can be proved by mimicking the proof of Theorem 14.51 given below for a larger 
class of domains X. However we give a proof of estimate (14.131) in Section El 

Note that the left-hand side \\g(k Q — k )\\ 1 of (14.121) can be replaced by ||fco — fcoll „-i,^ 

H 2 [Zj ' H 2 yZ) 

since g" 1 G C°°(Z) and the operator / — > f} -1 / is bounded in H~^{Z) for any open convex 
subset Z (with C 1 boundary) of X which satisfies Z C\ X. 

Under the assumptions of Theorem 14.41 and additional regularity assumptions on (a, k) one 
obtains stability estimates similar to those given in Corollary 14.61 given below for a larger class 
of domains X. 



4.3 Uniqueness and stability estimates for more general domains X 

Theorem 4.5. Assume that the open subset X ofW 1 is convex with a real analytic boundary 
and that inf^/y^r- S(x', v') > and mf^ XjV - )&r+ W(x,v) > 0. Let (cr,k) and (a,k) satisfy 
condition (I3.23p . Assume also that a, a do not depend on the velocity variable (a{x,v) = cr(x)) 
and k(x,v,v') = ko(x)g(x,v,v') and k(x,v,v') = k (x)g(x, v , v '), g(x,v,v') > 0, for (x,v,v f ) G 



X x 



5n— 1 



x § n where g is an a priori known real analytic function on X x S n x 



and 



where suppA: U supp/cn Q Z , (k , k ) G L°°(Z). Then estimate (14.111) still holds and there exists 
C 4 = C 4 (S, W, X, Y, Z, M, g) such that 



I An 



< CM 



+ 



~ 1 
til 2 

n W oo 



\k\\no\\A 



s,w 



A 



S,W\\r),T 



\z — z 



3 — n 

1 2 1 



Ti -fi)(r, z, z') 



(4.14) 



where M = max(||er| 



L°°(Y), \\ '\\l° o (Y)j 



Theorem 14. 51 is proved in Section El Assume that X is convex with a real analytic boundary 
and that inf( x /y) gr _ S(x',v') > and inf( a . )V ) e r+ W(x,v) > 0. Let Y and Z be open convex 
subsets of X, Z C X, Z C Y C X, with a C 1 boundary. Let g be an a priori known real 
analytic function on X x S"" 1 x S™" 1 , g(x, v, v') > for (x, v, v') G X x S™" 1 x S™" 1 . Let n > 0, 
r 2 > 0. Consider the class 

AT : = {(a, fc) G Ft+ r i(y) x L°°(Z x § 



n-l 



X 



0" 



<M l5 



k = k g, suppfc Q Z, 



\ k o\\ H %+r2 (z) <M 2 ). 



(4.15) 
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Note that there exist a function D\ : N x (0, +00) — > (0, +00) such that 

IM|l°°(y) < Di(n,ri)\\a\\ H n +n < D 1 (n,r 1 )M 1 , 
\\k \\L°°(z) < D 1 (n,r 2 )\\k \\ H ^+r 2 ^ z) < Dx{n,r 2 )M 2 , 
\\k\\L^{z) < \\g\\L°°(z)\\ko\\L°°(z) < Di{n,r 2 )M 2 \\g\\L^{z), 

for (a, k) G N. We also use the interpolation formula 



(4.16) 



H s (0) 



< 



a 2-"l II f || s 2- s l 

iPi(0)IU \\H s 2(0)i 
1 n 



(4.17) 



for si < s < s 2 and for (O, S\, s 2 ) G {(Y, — ~, | + n), (Z, — |, § + r 2 )}. Using Theorem 03] and 
f!4.16p . and applying (14.1 7p on / = a — a and / = k — k we obtain the following result. 

Corollary 4.6. Let (a, k), (a, k) G N. Then, for — ~ < s < ^ + 7*1 and /or < r < r 1; tnere 
exists C 5 = C 5 (S, W, X, Y, Mi, r 1 ,s) such that 



\W - (t\\ H s(y) < C 5 \\a 



<J lll°°(y)ll^s , ,vi/ — As.vkII^T' 



(4.18) 
(4.19) 



where («,«/) = and C 6 = CpL^r)^ (£>i(n,r) zs denned 6y flHB])j. 

In addition, there exists C 7 = C 7 (S, W, X, Y, Z, g, Mi, r 1; M 2 , r 2 , s) such that 



\ko — ko\\H"(z) — C*7 1| ^0 — ^0 
+ 



2 



A S w - A 



S,W\\t),T 



(4.20) 



\ko ~ ^o|l#f + r (z) — ^ 



A S w - A 



ir^\z-z'\) — (Ti-ri)(T,z,z') 

( t -\z-z'\)^{Ti-Ti){t,z,z') 



S,W\\r,,T 



for — I < s < I + r 2 and < r < r 2 , where (k, k!) = ^ 



n+2{r2 — s) 2(r'2— r) 
n+l+2r 2 ' n+l+2r 2 



(4.21) 
and (78 = 



C^DiinM 



(Di{n,r) is defined by (j4.16p ). 



Remark 4.7. (i.) Theorem 14.51 and Corollary 14.61 remain valid when: X is only assumed to be 
convex with C 2 boundary; the weight d D defined by (13.241) (resp. the function g which appears 
in the assumptions of Theorem 14.51 and Corollary 14. 6p is sufficiently close (in the C 2 norm) to 
an analytic weight #o,a on the vicinity of Z x S n_1 (resp. an analytic function g a on the vicinity 
ofZxS^x S"- 1 ); see proof of Theorem KB and [12l Theorem 2.3]. 
(ii.) When n = 3 then under hypothesis (13.231) . we have 



(r-\z-z'\)^(Ti-fi)(r,z,z') 



5> 

m=l 



m,S,W 



-A 



m,S,W, 



C(L l ((Q,7 j )xdX),L a °((0,T)xdX)) 



where the distributional kernel of the bounded operator ^^^(An^w/ - Aw,s,w) from -^ 1 ((0, n) x 
dX) to -^ 1 ((0, T) x dX) is given by r\ — r\. Therefore when n = 3 and under condition ( 13.231) . 
the right-hand side of the stability estimates (14.141) and (14.201) can be expressed with operator 
norms only (instead of using a norm on the distributional kernel of the multiple scattering). 
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5 Proof of Theorems I3.4L 13.5 



Proof of Theorem \3.5[ For the sake of simplicity and without loss of generality we assume 
v = (1,0, .. . ,0). Assume that condition (13.231) is satisfied. For n > 2 consider the following 
open subset of (0, +oo) x S™- 1 x S"" 1 



V := {(s,v,v') G (0,+oo) x S^ 1 x S£Ti | s G (0,r_(x,i;))}. 
Then we introduce the bounded function \l/ n on T> defined by 



(5.1) 



m-2 



W(x, f )(z/(x) ■ v)E(x, x — sv, x' )k(x — sv, v', v)S(x' , ?/)|z/(a; ) ■ v'\, (5.2) 



for (s,v,v') G T>. Note that from convexity of X it follows that t± is continuous on r T and 



fg a(x—pv,v)dp—f a(x—sv—p 



X — Xn—SV X — Xfx — SV 



E(x, x — sv, x'q) = e 
r_(x,v). Under (I3.23P we obtain that 



)d P 



\x sv\ \x Xq sv\ for V G 



jn— 1 



and < s < 



ty n (s, v, v') = for (s, v, v') G (0, +oo) x S%~+ x S£~_ such that x - sv £ Z, 

and the function \l/ n is continuous at any point (s, v, v') G V such that x — sv G Z. 

We first prove (EOoD for n = 2. Let r > t . From flPfl . (EUOD . it follows that 

1 



(5.3) 



Ji(t,x,x' ) 



TT-a 



/ r\\ <^(0,T_ (x,V )) \S ) ? ^ J 1? — (cos Q,sin f2) 

ao T -t COs(fi) J v , t n (l,0)-«, 



T — S 



' 2(r-iQ cos(fi)) 



= 7i,i( r > x > 4) + 7i,2(r, x, 4), 
where S™"^ 1 = {(cosf2, sin £7) | — a < Q < n — a } (0 < a < it) and 

X(0,7T-a )(^) 



7i,i{t,x,x ) 



7i )2 (r,x,x' ) 



We shall prove that 



o r - t cos(fi) J „,_ t n (i,o)-« 



T 2 -t 2 

in 



~ 2(t — tQ cos(fi)) 



_ n T - t COs(ii) L J p ,_ t n (i,o)- OT 



(5.4) 



(5.5) 



(5.6) 



- ' -'-o 

2(i — 1 () cos(f2)) 



y/r - toli,i{r, x, x' ) 

9bS T y tn 



W(x,v )S(x ,v )(v(x) ■ v )\u(x' ) ■ v \E(x' ,x) f to k(x - sv ,v ,v ) 



-ds, 



b 5 



(5.7) 



for i = 1,2. Then adding ( 15.7ft for i — 1 and z = 2, we obtain (13. 26ft . We only prove ( 15.71) 
for i — 1 since the proof for z = 2 is similar. Let r > to- Using the change of variables 

s 



- — ; tttt — 2L rr a , we obtain 

2(t— to cos(SZ)) 2 ' 



1 



7i,i(r,a;,ar ; 



*° , n , ,,X(o,r_(^(^)))( s + V a )*2(s,v(s,r),t; , (s,r)) 



o ■'o 



(5i 
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where 

--'-/ 



v(s,t) = (cosf2(s, t), sin ft(s, r)), f2(s,r) = arccos 
t„(l,0)- (s+2^)v(s,r) 



r 



2_ + 2 




2s+T-t 



t 







(5.9) 



i/(s,r) 



r+tg _ 
2 * 



Let s G (0,t )- From O, it follows that 

v(s, t) -> (1, 0) as r -> u'(s, r) -> (1, 0) as r -> (5.10) 

Note that using the definition of t> and using the assumption a; + e(x — x' Q ) G X for some 

e G (0,1) we obtain t = t_(x,v ). Note also that the function s t— > , 1 , s G (0,t ), is 

integrable in (0, to)- Therefore, using (15.31) . the boundedness of \l/ 2 on V and the Lebesgue 
dominated convergence theorem, we obtain (15 .7p . This proves (13.261) when n — 2. 
Let n > 3 and prove (15351) . From (E2D and (13TT0|) . it follows that 

7i(t,x,x' q ) = — ° ° ,„ 4 X(o,+oo)(y(?) ■ v)Vn(s,v,v') t n v n -sv dv, (5.11) 

J§n-1 tofO — Tf T s 

7" —in 



lo 

2(T-t vv ) 



for r > \x — x' Q \. 



Let $(fi, w) = (sin fi, cos(fi)a;i, . . . , cos(fi)c<j ra -i) for G (— | , |) and u = (ui, . . . , oo n -i) G 
gn-2 Tj s i n g spherical coordinates we obtain 



77 2 ( r-t sin(Q))" - 3 

— t 2 (*§ + r 2 -2t Q Tsm(Q)y 



Ti(r,x,4) = / cos(^)- 2 , 2 ^ (5.12) 



X(o,+oo)(^(a;) • $(Q,a;))* n (s, a;) , i/) tn«n-.»(n, W ) tftjrifl, 

s ~ 2(r-t sin(f2)) 



for t > tn. Performing the change of variables "r = o , T ~'° - on the first integral on 

u ° ° 2(t— to sm(SZj) 2 ° 

the right-hand side of (15.121) . we obtain 



= 2 2 - n tl- n (r 2 -tl)^ I ^:^Z 2( ^_ r)n _ 2 (5.13) 



t0 v/K^ - r)™ 3 



T'J V — 



fT 2_ t 2) dudr. 

n=arc S in(t- 1 (r T )) 

° 2(r+^21) 

„/. «n«n-»*("i") 



Therefore using (I5.13p . (15.31) and (15.21) and using Lebesgue dominated convergence theorem, we 
obtain (I3.29p . This concludes the proof of Theorem 13.51 □ 



Proof of Theorem \3.4\ For the sake of simplicity and without loss of generality we assume 
Vq — (1, 0, . . . , 0). Assume that condition (13.221) is satisfied. We consider the measurable function 
ty n defined by (15.21) for all (s,v,v') G V where V is defined by ( 15. ip . Under (13.221) we obtain 
that 

the function is continuous at any point (s, v, v') G V 

(i.e. for any (s, v, v') G (0, +oo) x S™" 1 x S™," 1 , x - sv G X), ^ 5 ' 14 ^ 
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The proof of (13.261) under condition (13.221) is actually similar to the proof of (13.261) under 
condition (I3.23p . Note that (I5.4p - (l5.6p still hold so that we have to prove that (15.71) still holds 
for i = 1,2. Again we only sketch the proof of (15.71) for 2 = 1. Note also that (I5.8I) - (I5.10I) 
still hold. Then using (I5.8I) -( 15TTU1) . (15 . 141) and (I5.2p and using Lebesgue dominated convergence 
theorem, we obtain (15.71) for i — 1. This proves (13.261) . 

Let n > 3. Formula (I5.12p still holds. Now assume that n = 3. We shall prove (13.271) . Let 
r > t . Using the change of variables «e = ^^^^^^m ^ the equa i ity <^> 
gives 

In (^±|) 

7i(r,x,4) = ^ ° (7i ; i(r,a:,x / ) + 7 1)2 (r, x, x' )) (5.15) 

where 

ji,i(r,x,x' ) = / / X(o,+oo)(y(x) ■ $(n(T,e),u))y 3 (s(T,e),$(Q(T,e),uj),v'(T,uj,e))dMde 
Jo Js 1 

(5.16) 




7i )2 (r,x,Xo) = / / X(p,+oo)(v(x) ■ $(£l(T,s),uj))y 3 (s(T,E) 1 $(n(T,e),uj),v'(T,uj,e))dujde 

(5.17) 



si 



and 



r 2 +t 2_ (r _ to) 2(l- £)(r + to)2£ 



Q(r,s) := arcsin^ * ^- v 7 j , (5.18) 

r 2 - t 2 

S(T,£) := 2(r-t 8in(fi(r, e )))' (5 ' 19) 

// x t ^o - s(r,e)$(n(r,e),a;) 

u(r,w,e) := ? r , (5.20) 

T-s{T,e) 

for < e < 1, u G S 1 . 

We shall give some properties of f2(r, e), s(r, e) and v'(r, u, e) for < e < 1 and a; G S n_1 . 
From (EHBD , it follows that 

fi(r, £)-»•-, as r ^ 4, for all £ G (0, 1). (5.21) 
From (EHHD and (QUI) , it follows that 

<r, e) = , r ^ + r) (r i <b) r ^ ^> £ ) = ST ( 5 - 22 ) 

V ^ (to + ^Cr-to^ + Cr + to) 25 V ^ t + r+(r-t o y-^(r + t )^ V ; 

for e G (0, 1). From (15T22D it follows that 
s(r,e) -> + as r -> t£, when £ G (1/2, 1), s(r,e) ^ t as r ^ tj, when e G (0, 1/2). (5.23) 
In addition, from (j51E) - (l5^Dj | . it follows that 

-r(r + t ) 2 (r - t ) 2£ + 2t 2 r(r + t ) 2e + r(r + t ) 1+2e (r - to) 



v'(t,u,e) 



2t r 2 (r + t ) 2£ 



v /_( r _ t() )4 £(r + t()) 2 _ (r _ to) 2 (r + to) 4 £ + 2(r 2 + t 2 )(r _ to)2£(r ~ to) 2e 

(r + to + ir-toY^ir + t^ 



x v ' u - r ~ ' uy ^ (5.24) 
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for < e < 1. Therefore 

v'(t, u, e)^v = (1, 0, 0) as r -> t+, (5.25) 

for < e < 1. 

Using (15.161) . (15.211) . ( 15. 23ft . ( 15.251) and Lebesgue dominated convergence theorem, we obtain 
7i i i(r,x,x[ ) ) lim^ 3 (s,v ,v )dujde 

r^t+ Jo Jgl s-+t Q 

= 2nE(x' ,x)W(x,Vo)S(x' ,v )(u(x' )v )(u(x)v )k(x' ,Vo,v ) (5.26) 

(we also used (15.21) ). 

Similarly, using ( 15.161) . ( 15.21ft . ( 15.23ft . ( 15.251) and Lebesgue dominated convergence theorem, 
we obtain 

j 1>2 {t,x,x' ) -> 2nE(x ,x)W(x,Vo)S(x ,v )(v(x )v )(v(x)vo)k(x,Vo, v ), asr^t + . (5.27) 
Statement E7J follows from fl57T5|> and HQS} - (lo~2Tj) . 

Let n > 4. We shall prove (13T28D . From (15TT2D it follows that 

7i(r,x,Xo) = 7i,i(r,x,x' ) + 7i j2 (r, x, x' ), (5.28) 

where 

' n 



2(t — £q sln(Q)) 

= /l ( » ) - (n) -' M^^S^ (5 ' 30) 

X / X(0,+oo)(^(^) • $(^, ^))X(0,r-( Z ,$(n,a/))))(s)^n(s, Ct>), f') , t n « n -^(n,^) dudQ. 

-"X 



2(r-t sin(S2)) 

First we study 7^1. Note that 

r - 1 mngji) r - 1 smgjj _v 

t§ + r 2 - 2t r sin^) tg + r 2 - 2t r sin(ft 2 ) ' 1 ' J 

for — I < < f2 2 < f and for r > t - Therefore using also the estimate cos(f2) < 1 we obtain 



(t 2 + r 2 -2t rsin(fi))«- 2 



- C0S ^J t 2 + r 2 - 2t r sin(fi) V*o + ^ " 2t r sin(fi) J " 2t r- 2 ' 1 J 

for n G (-f , f ) and sin(fi) < 7} (we used (jOIj) with u Q 1 n = Q and "fi 2 "= & and we used the 
estimate £5 + r 2 — 2toTsin(fi) > 2£ot(1 — sin(i7) )) , where 

sin 2 (y) cos 2 (Q) 

C := sup p-r = sup — — . (5.33) 

^g(o,2^) 1 - cos(y?) ne(-^,f) 1 ~ sm ( fi ) 

18 



Using (I5.29p . (15.321) . (15.21) and Lebesgue dominated convergence theorem, we obtain 
71,1(^,0/0) — > 2 2 - n tl~ n I" C ° s{Q)n ' / X (i/(a;) ■ $(Q,u;)) lim * ft (s, w), u )du;dn 
= 2 2 - n t - m / lim V n {s,v,v )dv 



il— n ../ ■■• ■ ' ■ ' ' ' 



t n E(x' , x)S(x' , v )(h , (x' ) ■ v ) / W(x, v )(jy(x) ■ v )k(x, v 0} v )dv . (5.34) 

-1 1 — V • Vq 



Now we shall study 7^2 defined by (I5.30p . Note that using the convexity of X we obtain 
/ _ x-sv-xq gn-i ^gngygj. v ^ g"- 1 anc i s ^ (0, t_(x, v )). Therefore using the change 

\x—sv—x \ X 0' _ ' 

2Tfo-(T 2 +tg)sin(f2) „ Q / ^ /a / \ dQ _ Q 2 -*q) 2 cosfflQ , r 



01 variables sm[\i ) - T 2 +t -z_ 2toTSin{n) , " fc ( 2 , 2 M cos l iZ J d n' - ( T 2 + ^_ 24or sin(Q')) 2 " 
obtain 

/ M f rcsin (^) , n „ w _ 2 (r-tosin(^)r- 3 
7l ' 2(T ' X ' Xo):= | f C ° S(fi) (tg + r 2 -2t rsin(^))"- 2 ^ 

w es«-2 X(o,r_(o:,0(n(T,n')^)))( s ( r ' ^'))*n(s(r, fi'), $(^( r > v') v >=0fQ' -^duodD! , 

i/(x)-*(n(T,n'),")>o 

-^(a:^)-*(n',-u)>0 

where 

n( T, O) = arcsin P7 = f + $ ^ ) , ,( T , ff) = Z^izj^Q, (5.36) 



r 2 + t 2 _ 2t rsin(fi / ) / ' 2(r - t sin(fi')) 

for fi' G (-f,f). Note that 

Q(t,Q') - |, $(0(t,0'),o;) - Uo , (5.37) 

s(r,fi') -> t , (5.38) 



for (Of, uj) G (-§, f ) x § n " 2 . Note also that from fl5T37D . it follows that at fixed v' = $(fi', -w) G 
§^,~i the condition X(o,T_(a:,$(n(T,fi')><«0))( s (' 7 ~> ^0) = 1 (which is equivalent to x— s(r, 0')$(f)(r, fi'), a>) G 
X due to convexity of X) for $(f2(r, O'), u;) G S™^ 1 is satisfied when r — to > is sufficiently 
small. Therefore from f!5.32j) (with "fi" replaced by "f2"') and from (15.351) and Lebesgue domi- 
nated convergence theorem we obtain 

7i, 2 (r,x,x' ) — ► 2 2 ~X~ n [ 2 r^O^ I X(o,+oo)(-^o) ' 
r^t+ J-TL i-sin^z; 

x lim ty n (s,v ,${Q',-u))dudQ' = 1 2 ~ n t\^ n I - lim ty n (s,v ,v')dv' 

s-^tg Jg*- 1 I — V-Vq s-yto 

= tk~ n E(x' , x)W(x, vq){v{x) ■ vq) [ S(x ,v')\u(x ) • v '| k(x, V ', v o)dv' . (5.39) 

Js n ~\ 1-V-Vq 

x a> 

Statement fl3~^8l follows from <^M) , flCTj) and (ICTjl . Theorem E31 is proved. □ 
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6 Proof of Theorems 14. 1L I4.2L 14.51 and Theorem 14.41 ( 14.131) 

Proof of Theorem \4.fy We now prove (14.51) . Let x' Q G dX. For e = (£1,62) £ (0, +00) 2 and 
£3 e (0,+oo) let (/ ei ,<fe a ) G C^dX) x C^M) satisfy 

9e 2 > 0,f ei > 0, suppfe C (0,min(e 2 ,»7)), (6.1) 
supp/ £l C {x 1 G dX I \x' - x' \ < £1}, (6.2) 



n 



g e2 (f)d1? = l, / / £1 (x')^(V) = l, (6.3) 



dX 



'(0,T)xdXx(0,r])xdX 

for £ = (£1,62) (0, +00) and £3 G (0, +00), where 



for £ = (£1,62) G (0, +00) 2 . Therefore £ := g £2 f ei is an approximation of the delta function at 
(0,x' ) G M x dX for £ := (£ 1; £ 2 ) G (0, +00) 2 . Let G L°°((0,T) x dX) be defined by 

^ e3 (t,x) = X (-ss,e 3 )(t- \x-x' \)(2 X(0 , +oo) ((E-E)(x,x' )) - 1), (t,z) G (0,T) x (6.4) 

for £ 3 > 0. From flUlD and fl3~T2l it follows that 

/ ifj £a (t,x)(A s , w ~ A s ,w)<Pe{t,x)dtdn(x) = I (if> ea , <f> e ) 

J(0,T)xdX 

ip £3 (t, x)(p e {t', x')(T% — Ti)(t — t', x, x')dtdf/,(x)dt'dn(x'), (6.5) 



r / , , x f , / x , / /, /x -fife a/) — E(x, x 1 ) 

Io(Aa,4e) = / ^ £3 {t,X)4) £ {t-\x-x'\,x') K \ ' „ \ ' 

x [W(x, f )S(x', f )(^(x) •f)|z/(x / ) • f |] _ a-a/ dtdn(x)diJ,(x') . (6.6) 

I x — a/ I 

From fl3~TBl) . (|5T7L (I530j) and fl3~2TD it follows that 

(r- |x -x'|) — (ri - ri)(r,ar,a/) G L°°((0,T) x^Ix dX). (6.7) 

Combining (ESI) and the equality ||0 £ ||li((o )7? )x9jo = 1 and the estimate \\ip £3 \\L^((o,T)xax) < 1 
and (16.71) we obtain 

Ws 3 , <Pe) < \\A s ,w - A s ,w\\v,T + CAi(Vt 3 , £ ), (6.8) 

for £ = (£i,£ 2 ) G (0, +00) and £ 3 G (0, +00), where C = ||(r — |x — a/)) - ^^! — ri)(r, £, x') 
\\L^((o,T)xdx x xdx xl ) and 

/ Ti — 3 

Ajfaks,^) = ((LT)txgXxX(0tV)tixgx , i> e3 (t,x)<t> e (lf,x')(t-t' - \x-x'\)—dtdfi(x)dt'dn(x). (6.9) 

|x-a:'|<t-t' 

Note that the function $ 1)£3 : [0, 77) x <9X — ► K. defined by 

/ n — 3 

$ l£3 (t',a/) := / ip £3 (t,x){t-t' - \x-x'\)~dtdfx(x), (t',x) G [0, 77) x dX, (6.10) 

J (0,T) t xdX x 
\x— x'\<t— t' 

is continuous on [0,i]) x dX for £ 3 G (0, +00). Therefore from fl6.1l) - fl6.3p and the equality 

Ai(^ £3 ,0 £ ) = f {Q>v)xdX (f) E (t',x l )®i !E3 (t',x')dt?diJ l (x') it follows that 

lim lim lim Ai(^ £3 ,0 £ )= lim $ lj£3 (0,a/ ) = (6.11) 

£3— >0+ £2^0+ £1— >0+ £3^0+ 
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(we also used (16.101) . (16.41) and the Lebesgue dominated convergence theorem to prove that 
lim £3 ^ + $x )£3 (0, x'q) = 0). Note that under condition (13.231) the function $o,£ 2 ,£3 : — > M. 
defined by 



(0,T) t x8X x 
x — x'\<t 



ipe 3 {t,x)g £2 {t- \x-x'\) 



E(x, x') — E(x, x') 



/ln-1 



x [W(:c, v )S(x', v)(u(x) ■ v)\v(x) ■ v\] v _ x-x> dtd/i(x) 



is continuous on dX, for (e 2 , £3) G (0, +00) 2 . Therefore from the equality Io(i/j e3 , 4>e) 
X f ei (x')dfi(x') (see (16.61) ) it follows that 



(6.12) 

lax ®o,s 2 ,e 3 {x') 



lim I (ip t 



£1^0+ 



£35 "re) 



$o, £2 , £3 ( x o) ; for (^,£3) e (0, +00) 



(6.13) 



Therefore using the Lebesgue dominated convergence theorem and (I6.12p we obtain 



lim lim lim I (ip £3 ,<f) £ ) 

£3^0+ £2^0+ £1^0+ 



E(x, x'q) - E(x, x' ] 



dX x 



\x — x: 



1 ln-1 



[W(a;, v )S(x' , v)(u(x) ■ v)\u(x' Q ) ■ v\ 



dfi(x). 



(6.14) 

Combining ( 16. 14ft . ( 16.111) and (16. 8ft we obtain the formula (14.11) . Using ( 14. ip and the estimates 
inf (a;>')er_ S(x', v') > and inf( x ,v)er + W(x,v) > and the change of variables x = x' Q + 
t + (x' ,Vq)vo ( i _ i?n-i dii{x) = dvo) we obtain ( 14. 5 h where the constant C\ which appears on the 

right-hand side of (14.51) is given by C\ = [ini( x ',v')er- S(x',v') inf( Xjl) ) gr+ W(x, v)) 

We now prove (14.61) . Let x G dX be such that px' + (1 — p)x G Z for some p G (0, 1). We 
set t = \x - x'\ and v Q = ^4. From (l3~T6l) . (l3~T71) . (l3~2UD and (EOT!) it follows that 

(r- |x - x'oD^^i -^K^x^q) < (r- |x - a^D^ra - f 2 |(r, z, 2;') (6.15) 

3 — n ~ 

+||(s -\z- z'\) — (T 1 - T 1 )(s,2;,z / )IU-((o,T) s xax z xax z ,) 5 

3_ n 

for r > |x — x' |. From ( 13. 17ft . (1 3 . 2 1) ( T3 . 2 1 1) it turns out that hm T ^ x _ x > Q \+(T — \x — x / |) _ 2 _ |r 2 — 
r 2 |(r, z, z') = 0. Therefore applying (13.261) and (13.291) on the left-hand side of (16.151) we obtain 



2 2 n \x-x' \ n 2 C n S(x' ,v Q )W(x,v )\v(x ) ■ v \(u(x) ■ v ) 



X 



P 2 (*0~P) 2 

3— 71 ■— 

< \\(S -\Z- Z'\) — {V X - T 1 )(s, Z, Z?))\\L°°«0,T) a xdX z xdX z ,), 



dp 



(6.16) 



where C n = 2 if n = 2 and C n = Vol„_ 2 (§ n 2 ) if n > 3. Then note that Cx '■= infmeex, zez v { x i)' 
> since X is a bounded convex subset of lR n with C 1 boundary and Z G X. Therefore 



( 14.61) follows from ( 16.161) where the constant C 2 which appears on the right-hand side of ( 14.61) 



is given by C 2 



n— 1 n— 1 

2^~diam(X)^~ 



C n C\ ini {x ,y )GT _ S(x',v') inf( a . 1) „ l)er+ 



. Theorem 14.21 is proved. 



□ 



Proof of Theorem \4.1\ We prove (I4.ip . Let x' G <9X. For e = (e 1 ,e 2 ) G (0, +00) 2 and e 3 G 
(0,+oo) let (f £1 ,g £2 ) G C^dX) x C\R) satisfy flSHJ-Q and tp £3 be defined by (Q. First 
note that ( l6.5l) -( l6T6i) still hold. From Theorems 13.21 and 13.31 it follows that 



\x-x'\ n -i(T- |z-a/|)*(ri-fi)(T,a;,z') G L°°((0,T) x dX x dX). 



(6.17) 
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Combining (16. 5p . (16.171) and the equality H^eIIli^o^xm:) and the estimate ||V'£3lU oo ((o,T)x ( 9x) < 1 
we obtain 

We 3 , <t>e) < \\A S ,w - A S ,w\\v,T + CA 2^e 3 , <t>s), (6-18) 

7 3 ~ 

for e = (£1,62) G (0, +00), £3 G (0, +00), where C = \\\x — x'\ n ^i(r — \x — x'\)^(Ti — r 1 )(r,x,x / ) 
\\L°°((a,T)xdx x xdX x , and 



A 2 (V £3 ,0 e )=/ fj^x^if^) dtd^{x)dt'd^{x'). (6.19) 

|i-i'|<i-i' 11/ 

Note that the function $ li£3 : [0,T) x dX — > R defined by 
<W*V):= / 9 , 1 „n 1 -^jdtdfi(x), (t',x')e[0, V )xdX, (6.20) 

J(°,Thxdx x \ x _ x!\n-z(t - t 1 - \x - X'\)* 



\x-x'\<t-t' 



is continuous on [0, rf)xdX. Therefore using (16.191) . (I6.1l) - (l6.3p and using the equality A 2 (ip £3 , 4>e) 
f(o, v ) t ,xdx x , ^3,e 3 (t' , x')(j) £ (t' , x')dt'd/j,(x') we obtain 



lim lim lim A 2 (^ 3 ,0 £ ) = lim $3, B3 (0,a/ ) = 0. (6.21) 

£3^0+ £2^0+ £1— >0+ £3— >0+ 

(we also used ( 16.20(1 . (16.41) and Lebesgue dominated convergence theorem to prove lim e3 ^ + ^3,£ 3 (0, x' Q ) 
0). Note that under condition ( 13.22(1 the function $4,£ 3 ,£ 2 : dX — > R defined by 



a> , f i (± \ (± 1 n\ E dXxdx{x, X 1 ) - E dX xdx(x,x') 
$ 4,£ 3 ,£ 2 (x) = / ^£ 3 (t,Z )& 2 (t- \X-X\) — — 

J (0,T) t x3X x I >• — >•' I ' 1 



\x — X 



\m—x'\<t 

x[W(x,v)S(x',v)(v(x)-v)\is(x')-v\} v=x - x , dtdfi{x), (6.22) 

I x — x 1 1 

is continuous on <9X. Therefore from ( 16. ID — (16.31) and the equality Iq{i^ £3 ^ 4>e) — fg X fa { x ')$ 'i,e 3 ,e 2 ( x ')dn>(x') 
(see ( 16.61) ) it follows that 

lim Jo(^£ 3 ,0e) = $4,£3,£2K)> for ( £ 2,£s) e (0,+oo). (6.23) 
£1^0+ 

Then using Lebesgue dominated convergence theorem and (16.221) we obtain 

r E 9X xdx(x, x' ) — E 9X xdx(x, x' ) 
lim lim lim / o (V £3 ,0 e ) = / , j— j (6.24) 

£3^0+ £2^0+ £1^0+ Jq Xx \X-X' \ nl 

x [W(x, v)S(x' ,v)(u(x) ■v)\u(x' ) -v\] x - x t dfi(x). 

v= ; _ ¥-. 

I X X Q I 

Combining (|E23j) . (I6T2T]) and (I6~TS|) we obtain (fl~Tj) . 

We prove (I4.2p -( l4~4"l) . Let x G dX be such that px' + (1 — p)x G X for some p G (0, 1). Let 
/3 n : {(r, 2, z') G (0, T) x <9X x dX \ z ^ z', r > \z - z'\} -> R be defined by 



r 2 — |z — z'| 2 , if n = 2, 



/3 n (r, z, z') = I if " = 3, (6.25) 



ln( t 

r|z - ^T" 2 , if n > 4, 



for (r,^, z') G (0,T) x <9X x dX, z ^ z', r > |z - z'|. From fl3TT3l) - fl3TT5l) and (l3TT71) -( l3TT9l) it 
follows that 

/3 n (r,x,x / )|(7i - 7i)(r,x,Xo)| < (3 n (r, x, x' )\(T 2 - T 2 )(t, z, z')\ (6.26) 

+ HA1O, 2, z')( r i - f i)(s, 2, 2 , )IU ao ((o,T).x8Jr,x9x^)- 
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From (16.251) and (I3.17I) - (I3.19I) it turns out that \im T ^\ x _ x > \+ /3 n (r, x, Xo)|(r 2 — r 2 )(r, z, z')\ = 0. 
Therefore applying (l3~26l) (resp. fl3~27j) . (JS2HD) on the left-hand side of (IB^o]) we obtain (TO) 
(resp. (14.31) . (14.41) ). Theorem 14.11 is proved. □ 

Proof of Theorem \4-5[ We first prove (14. lip . We extend a and a by outside Y. For a bounded 
and continuous function /on7 consider the X-ray transform Pf : § n_1 xK n ->K defined by 
( 14.71) (we extend / by outside Y). We recall the following estimate 



<( f [ \Pf(v,x)\ 2 dxdvY , (6.27) 



jn— 1 



e -Jo o-(z +ot,i;)<2s _ g - J T u <r(iE +s'u,r>)ds 



where IT, := {x G W 1 \ v ■ x = 0} for v G S n . Note that using the estimate \\ct\\oo < we 
obtain 

rT + (x' Q ,v) 

/ a(x[, + sv,v)ds < Mt + {x' ,v) < Mdiam(X), for (x' ,v) G T_. (6.28) 

Jo 

Replacing a by a on the left-hand side of (16.281) we obtain an estimate similar to (16.281) for a. 
Therefore using the estimate \e h -e* 2 | > e" Mdiam(x) |ti -t 2 | for t 2 ) G [0, +oo) 2 , max(*i,t 2 ) < 
Mdiam(X), we obtain 

> e- Mdiam W\P(a-a)(v,x' )\, (6.29) 

for (x' ,v) G r_. Integrating the left-hand side of (14.51) over dX and using (I6.29p . we obtain 

J \P{o - a)(v,x' )\ d£(v,x' ) < e Mdiam ^Vol(dX)C 1 \\A s ,w - As,w\U,t, (6.30) 

where C\ is the constant that appears on the right-hand side of (14. 5p . Note that using that X 
is a convex open subset of W 1 with C 1 boundary we obtain j T \P(o~ — a)(v,x' ) \ d£(v,x' ) = 

fgn-i Jn \P( a ~ cr)(v,x)\dxdv. Therefore using (I6.30p and the estimate \P(cr — a)(v,x)\ 2 < 
|| cr — cr|| i oo(y)diam(X)|P(cr — a)(v,x)\ for (v,x) G TS n_1 (see (16. 28ft and the estimates a > 0, 
cr > 0) we obtain 



i 

2 1 1 



P(cr - a)(v, x)\ dxdv ) < C 3 \\cr — cr\\£o\\A StW — A SiW \\* T . (6.31) 



where C 3 = (diam(X)e Mdiam ( x )Vol(9X)d) 2 . Combining (IOTP and (IOTP we obtain pi) . 

We now prove (14. 141) . Let / G L 2 (X), supp/ C Z. We consider the weighted X-ray trans- 
form of /, P#f, defined by 

P$f(x,v) = I f(pv + x)$(pv + x, v)dp, for a.e. (x, v) G r_, (6.32) 
Jo 

where d : X x S n_1 — > (0, +oo) is the analytic function given by 

•&{x,v) = (T^(x,v)r + (x,v))~ IL 2 1 g(x,v,v), for (x,v) G X x S n_1 . (6.33) 
From [T2|, theorem 2.2] and from [T6], theorem 4] we obtain 

H - i(z) <C\\Pof\\ L2{r _ M) , (6.34) 
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where C = C(X, Z,g) is a constant that does not depend on /. Let x' G dX and let x G OX 
such that px' + (1 — p)x G Z for some p G (0, 1) where Vq = ^.J^\ an d = \x — x' Q \. Note that 
using (13.231) (since k G L°°(Z) and suppA; C Z C {x G X | inf^ax |x — x'| > 5}), we obtain 



p 2 (r+(4,^) -p) 2 



^rdp < \\k\\ L ™ {Z ) 



t+(x ,v' )-6 



p 2 (t+KX) -p) 



< 



wS-^T+tfoA) < ||fc||L-(z ) r^ 1 )diam(X). 



We use the estimate 

\P^k -k )(x' ,v' )\ < e p ^'o^\Pik (x' ,v' )\ 



-Pa{v' ,x' ) _ p -Pa{v' ,x' ) 



EI d P 

2 

(6.35) 
(6.36) 



+e 



Pa(v' Q ,x' Q ) 



- p ^'o^PMx' ,v' Q ) - e- pa ^Pik (x' ,v' ) 



Integrating both sides of inequality (I6.36P over v' G $ n , * and using the estimate e Pa ( v '°' x °> < 
e Mdiam(x) ) and uging ([^d (|^5]) - (|Qj) . W e obtain 



P*(fc - fc )IK,^)|i/K) • v \ dv < 5-("- 1 )diam(X)e Mdiam WC 1 ||fc|| 00 ||A 5iW - A SiW ||,, T 
Vol(§ n - 1 )e Mdiam ( x )C 2 



n-1 
I 

0' 



(r-lz-^D-rCTi-rOCr.z,^) 



L°°((0,T)xdXxdX) 



(6.37) 



where Ci and C% are the constants that appear on the right-hand side of (14.51) and (I4.6p . 

From the estimate \P$(k — k )(v' Q , x' )\ < \\k — fcllxoo^^^^diampT) for a.e. (x' ,v' ) G T_ 
(see (I6.35P ). it follows that 



\\P«(ko-ko)\\h<r ^ ) <\\k-~k\\ L ^( Z )d- in - 1) dia l m(X) [ [ \P4k -~k )(x' ,v' )\\u(x' )-v\dvdfi(x' ) 

JdX Js^ 1 

x o> 

Combining flQ7l) -f l6T38l) and fOIft we obtain fT4~T4l) . 



(6.38) 
□ 



Proo/ o/ Theorem gg(B3SD. We first prove fl6\42l) given below. Note that from fl3T24D -f l3T25l) . 
it follows that 

-T+{x,v) 



\P*of(v,B)\ < 



1 



t(r+(x, v) — t) 



dt = C\\f\\ t 



(6.39) 



-dt. 



for (v, x) G T_ and for / G C(X) n L°°(X), where C = £ -^L- . 

Note that u(x) = x and z/(x) • (a; — x' ) = \v(x' Q ) ■ (x — x' Q ) | for (x, x' ) G dX = S 1 . Therefore 
from (jO|) . it follows that 



< c 


\/ T 2 — : - :' 







W( x 'o)- v 'o\ 2 E(x,x' )P A) k (v' ,x) - E(x,x' )P# Q k' (v' ,x) 

(6.40)' 

for (x,x' ) G <9X 2 , x ^ x' and ^ = (we also used P# o k Vo (v ,x) = g(vo,v )P^ k (v , x) 

and the similar identity for k), where C = 2 i n f r Smf r Wmi Tgivv)' ^ n addition, from (16.401) . 

(16.36P (with P$ in place of "P?"), and from (I6.28|) and (I6.39P (with k in place of /), it follows 
that 



|^' )-^| 2 |P^(^o-^o)K,^)| < e Mdi ^C\\k \ 



e -Pa(v' ,x' ) _ e -Pa(v> ,x> ) 



+e 



Mdiam(X) 



\p(x' )-Vq\ 

(6.41) 
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for (x' , v' ) G T_ (we also used the estimate \v(x' ) -v' \ < 1). Performing the change of variables 
u x = x' + t + (x' , v'^Vq" ( i^^n^rr dpi(x) = dv'o) on the left-hand side of ( 14. ip . we obtain that the 
estimate (T4~5D still holds. Using ([43]) . (ET4T]) and (15~4T?]) (and (T4~%]) ). we obtain that there exists 
a constant C" such that 



n-l 
I 




z/(a/ ) • Uq| I P(p(A;o - A; )) (rf 0) v' )\di/ < C" (||fco||°o \\A s ,w - A s , w\\ v 



V 2 - |z — -2'| 2 (ri — Ti)(t, Z, z) 



(6.42) 



for x ; G dX. Moreover, using (I6.39j) (with k — k in place of "/") and Cauchy-Bunyakovski- 
Schwarz estimate, we obtain 



\P(p(k - ko))\ 2 (v ,a/ )d£(xk,v ) 



<C\\k -k 



Olloo 




dX JS 1 , 



\P(p(k - k ))\(v , x' )\v Q ■ u(x' )\ dv' dp(x' ) 



<Cpo-^o||tV27r f / / \P(p(k -k 

\ JdX is 1 , 



o))|(vo,^ )l v o • K^l^odM^o) • (6-43) 



Finally combining (16T42|) - (16T13| ). (IQTj) (and the identity f r P(p(k - k )(v,x' o ) d£(v,x' ) = 
f s „-i J u \P(k — k )(v , x)\ 2 dxdv), we obtain (14.131) . □ 



7 Proof of Theorem [372 

For < 6 < a we remind that 



2- 



2tt 



(7.1) 



, a-6sin(fi) Va 2 - b 2 

We will use the following Lemma 17.11 to prove Theorem 13.21 ( 13.131) , ( I3.14p and ( I3.15P . 

Lemma 7.1. Let n > 2. Let N denote the nonnegative measurable function from (0, T) x dX x 
R n to [0,+oo[ defined by 

at/ A /i /ix f (r - (x - x') ■ v)"~ 3 , 
N{t,x,x ) = X(o,+oo) [r — \x — x \ ) I — : — , dv, 



-i \X — X — TV 



2n-4 



(7.2) 



for (r, x, x') G (0, T) x dX x W 1 . When n = 2, then 

N(r,x,x') = X(o,+oc)(r -\x- x'\ 

for (r, x, x') G (0, T)xdX xW 1 . When n = 3, then 
N(t,x,x') = 2tt 



2tt 



T* — \X — X 



/|2 ' 



X(0,+oo)(r - 


X 


-x'\) 


T 


x — x' 





111 



T — \x — X'\ 



for (r, x, x') G (0, T) x dX x M n . Wften n > 4, £/ien 

r|x — a/| n ~ 2 iV(T, x, a/) < oo. 



sup 

(r,i,i')£(0,T)x3Xxl» 



(7.3) 



(7.4) 



(7.5) 
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Proof of Lemma\7l\ Let (r,x,x') G (0, T) x dX x R n . We first prove (Q. Let n = 2. Note 
that 



2tt 



N(r,x,x') =X(o,+oo)(r -\x - x'\) / . — — — dO. 

J t — \x — x I sm(iz) 

Therefore using ( 17. ip we obtain ( 17. 3D . We prove ( 17.41) . Let n = 3. Note that 

iV(r, x, x') = 2 ^(oW r -|z -^l) P d ^ , 2 _ 2 _ _ . ^ 

2r|x — ar | J_* all x ' 

which gives (17.41) . 

We prove (1731) . Let n > 4 and let (r, G (0,T) x OX x R n be such that r > \x — x'\ 
(we remind that N(t,x,x') = if r < |sc — a;'|). Using spherical coordinates, we obtain 



n-2\ I 2 V ~~ |J> — J> | amyaojj „„„/n\n-2. 



(r- 




sin(0)) n - 3 


( \x — x' 


2 + r 2 - 2r|x - x' 


sin(0)) n - 2 



iV(r,x,x / )=Vol n _ 2 (§ n - 2 ) 

Performing the change of variables "r = 2 ( r -|^"-a;'|sinffi)) ~~ zH£z£_l" ; we obtain 

x = Vol ra - 2 (S"- 2 )(r 2 -|x-a:f)^ r\ x ~ x '\ ^ r (\x - x'\ - rf' 3 ^ 





\X — X'\ n 2 Jq 


2Vol n _ 2 ( 


§n 2)( r 2_ Iz-z'^V /- 1 




x — x' n ~ 2 Jq 


2Vol n _ 2 ( 


§n 2 )|( r 2 _ | x _ x '|2)" 2 3 /• 




1 />"■ / > ■ / 1 ?t< 2 / 


2Vol n _ 2 


;§"- 2 )(t 2 - (x-x'H^r 1 



r-n— 3 



\/r(|x — x'| — rj 

ctr 



' T-\X-X'\ _|_ r '| rl - 2 j- T+|x-X / | 



r 



in-2 



1 



< ^ ir / : dr 



\x — x f \ 

" k-x'|™- 2 2 y (l^z£l + r )^± r 



n — 3 
2 



< ^Vol n _ 2 (§"- 2 )|* - x'| 2 -" f T+ ' X ^ 
n — 3 \ 2r 

on— 1 

< -Vol n _ 2 (§"- 2 )|a;-xf-V- 1 , (7.7) 

which proves (17. 5p . □ 

We are ready to prove Theorem 13.21 First we give an estimate on the simple scattering 
term. From (13. lOj) it follows that 

Mt,x,x')\ < 2 n - 2 \\W\\ 00 \\S\U\k\\ 00 I 1 (T,x,x') (7.8) 

for a.e. (r, x, x') G M x dX x dX, where 

h(j,X,x') = X(0,+oo)(t - \x - x'\) / Xsuppfc(^ - SV) _ T 2_\ x _ x/] 2 

J§n-1 l S 2(t-v(x-x')) 

(t-(x- x') ■ v) n ~ 3 

dv. (7.9) 



\x — X' — TV 



2n-4 
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Let (r, x, x') G (0, T) x dX x dX be such that x 7^ x' and r > |x — x'\. Assume without loss of 
generality x' — x = \x' — x\(l, . . . 0). 

First we prove (!333l - (l3TT5l) . From (J72D and (Q, it follows that 

Ii(t,x,x') <]V(t,x,x'). ( 7 - 10 ) 

Combining (|73D (respectively flUD, dUSD) with (J72D and (17TTUD . we obtain (1X13]) (respectively 

©H, (EUD). 

Now assume that k G L°°(X x S"" 1 x S"" 1 ) and 

supp/c C {x G X I inf |y — x\ > 5} for some < 5 < 00. (7-11) 

y£dx 



2 I / 1 2 

Let f G S n_1 and s := T ~r ~ g 1 . Straightforward computations give s + \x — x' — 
Using ( 17.1 ip we obtain that 



SV\ = T. 



if r < 5 or s > r — 5, then x — sv G" suppfc. (7-12) 
Using flUD and (17TT21) . we obtain 

if r< 5 then Ii(t,x,x') = 0. (7.13) 

We prove (I3.16P for n = 2. Using ( I7.10p . we obtain that Mr, x,x') < , 2?r = for r > <5. 

v<5i/t— |:r— x'| 

Combining ( 17.81) with this latter estimate and ( 17.131) . we obtain ( 13.161) for n = 2. 

We now prove ( 13.161) for n > 3. Let n > 3 and r > <5 (the case r < 8 is already considered in 
( I7.13P ). Performing the change of variables "r = 2 ( r -|3^a/|sin(n)) ~ T ^ X 2 X - " with "t> = $(f2, u;) : = 
(siii(n),cos(n)w), G (-f , f), uj G § n ~ 2 " on the right-hand side of (ESI), we obtain 



/iCr.x.ar') = 2 



2' 2- 

n— 3 

, t*- I .. >•' I - ^ — r- /" ' ./,./!,._ _ ,.\ 

2-n > 



(r 2 - 






x' 


) 2 






— x' 


n- 


-2 



\/r(\x — x'\ — r) 



■ r-\x-x'\ _|_ r ^ n -2 ^+\x-x'\ _ r ^ n ~2 



[XsuppfcOr - sv)] j^^i) dudr. (7.14) 



2 

Now assume r > | + |x — x'|. Then 
Is - x'?- f lX ~ X ' 1 Vr(\x-x>\-rr 3 dr 

J + r )n-2(I^£z£l _ r )„-2 

< ( — ] |x — x'| 2 ~ n / y/r(Jx — x'\ — r) dr = I - ) / a/ r(l — r)<ir < S 



Therefore using (I7.14p we obtain 



{r-\x-x'\)-^h{T,x,x') < 2"- 2 Vol ri ._ 2 (§"- 2 )(T + diam(X))^ ( v - ) . (7.15) 



4-2n 



Finally assume 5<r<2 + |^ — #'| and |x — x'| < r < T. From ( 17.121) . it follows that 

, , _ , Vol„_ 2 (§"- 2 )(T + diam(X))^ " + J'* ) ^(Ix - x'| - r^" 
(r— |x— x |) 2 I x [r,x,x)< ' 



2™- 2 |x — X'\ n ~ 2 J ( t-\x-x'\ 1 r -\ n -2C I±]£z^l r ')n-2 

(7.16) 
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where / 

r-[j,x,x):= , r + (r,x,x):= . (7.17) 

Note that 

r +( T ' x ' x '} \/r(\x - x'\ - r) n 3 r^r 1 \/r(\x - x'\ - r) n 3 

df — 2 / i 7] \ i 7] df 



,*,*') + r )n-2(I±^l _ r) n-2 J r _ {T!X>xl) (I=k=*l + r )«-2(I±l^l _ r )n-2 



r_(r,x,x') ^ 2 rTJ (_ 2 — '"J Jr-{r,x,x') y 2 1 ' / V 2 

|x— Jc I \x — x r I 

<2(-) \x-x'\ n - 3 / ; , cir = 2 n - 1 r 2 - n |a;-x , r- 3 / , T , dr. 

(7.18) 

Using (17.1 Tl) we obtain 

['^r 1 1 f In (f), if it = 3, 

pEEH + r) ^* = C <"' T ) : = { A ((I)-" - (5) 3 -) otherwise. <" 9 > 

From (J71SD, (EHHD, (173$]) and the estimates 5 < r < f + \x - x'\, it follows that 

(r-\x-x'\)-^h(T,x,x') < 2"Vol ri _ 2 (§"- 2 )r (n - 1) (T + diam(X))^C(n,T), (7.20) 

where the constant C{n,T) is defined in (17391) . Combining (ESD with (17331) . (17351) and (TT20D . 
we obtain ( 13.161) for n > 3. □ 



8 Proof of Theorem 13.3 



We shall use the following Lemmas 18.11 18.21 I8~3l and l8. 41 Lemmas 18 . H 18. 21 18311 and 18.41 are proved 
in Section [HI 

We introduce some notation first. Let m > 1 and z' , z e R n such that z ^ z' . Let fi > 0. 
We denote by £ m ,n(/i, z, z') the subset of (R n ) m defined by 

E m>n (fi,z,z') = {(y 1 ,...,y m ) e (K") m | \yi\ + ... + \y m \ + \z-z' -y x - ... -y m \ < //}. (8.1) 

When /i < |z — then £ m>n (fJ,, z, z') = 0. 

Lemma 8.1. Let J 2 be the junction from (0,T) x dX x R™ defined by 

J 2 (fi,z,z')= I r L^N( f Ji-\y\,z,z' + y)dy, (8.2) 

j £i in (ji,z,z') \y\ 

where N is defined by ( 17.21) . Then the following statements are valid: 

sup J 2 (//, z, z') < oo, when n = 2; (8.3) 

(M,z,z')e(0,T)x8XxR" 
^i>|z — z'\ 

sup (/i— |z— ^l) -1 //!^— z'\ (l + In | — — j j j J 2 ([i, z, z 1 ) < oo, whenn = 3; 

(0,T)xSXxR™ \ VA 4 ~~ |Z — £ | / / 

(8.4) 

sup (yit — |z — z'|) Vl- 2 ~~ z'| n ^O", z, z') < oo, when n > 4. (8.5) 

(M,z.z')e(0,T)x9XxR™ 
M>|z-z' 



(n,z,z')e(0,T)x8XxR" 
fj,>\z — z' 
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Lemma 8.2. Let m > 3 and let J m be the function from {(r, x, x') G (0, T) x W 1 x M. n | < 
\x — x'\ < t} to M defined by 

r ( i\ f dy 2 ... dy m -i 

J m (T,X, X) — / I : " : : ~, : -, (8.0) 



m — 2.n \ 



for (r, x, x') G (0, T) x K" x W 1 , < \x — x'\ < r. When n = 3, then there exists a constant C 
which does not depend on m such that 

\x — x'\ \ \T—\x — x'\JJ [m — 3J! 

/or (r, x, x') G (0, T) x 1" x M n ; < \x — x'\ < r . When n > A, then there exists a constant C 
which does not depend on m such that 

Jm(r,x,x') < CiT-lx-x'Wx-x'l 2 -"—^- 1 ^^ 



(m-3)! 



for (r, x, x') G (0, T) x lR n x R n , < \x - x'\ < r. 



Lemma 8.3. Let n > 2. Let (r, x, x') G K x IR n x R n be such that r > \x — x'\ > 0, the following 
estimate is valid: 

VtM£ljT , x , XJ) < vo...- 2 (s->(r + i,-,i) u^w y\ (89) 

where £i >n is defined by (18.1 1) . 

Lemma 8.4. Let £? fte £/ie function from {(/i, z, 2;') G (0, T) x <3X x 1R 3 | > \z — z'\ > 0} to 
M. defined by 

B(P, z, ;<):=/ In f"-|"| + l-— "'-■'I) * = / l n (»-M + \M»)-y\) 

Je 13 (n,z,z>) \l*-\v\-\z- z -y\J Je 13 Qt,u>{i,ofl)fl) \y\ - l(*o,0,0) -y\J 

(8-10) 

for (fi,z,z r ) G (0,T) x OX x W 1 where t = \z — z'\, z' ' , \i>\z — z'\. Then we have: 

sup 0- Iz-z 7 !)" 1 fl + ln ^"^h l ^(/i, z, z') < 00. (8.11) 



(M,z,z')e(0,T)x8Xx 
H>|z-z'|>0 



We also need the explicit expression of j m , m > 2, to prove Theorem 13.31 

j 2 (t,x,x'):= / / (u(x) ■ v)W(x, v) [E(x, x - (r - \y\ - s^v, x' + y, x') 

J vse l n (-i-,x,x ') Jgn-i 

x' + y£X x < + 



X(o,t_ («,«)) (t - 1 2/ 1 - si)/c(a; - (r - s x - |y|)t>, v u v)k(x' + y,v\vi)S(x\v') 
/a „i 2 n - 2 (r - lyl - (x - x' - y) • v) n ~ 3 , , 

a:' — y— (t — a-] — I y I )-» 



"1 = 



: 

/_ y 
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and 



j m {r,x,x'):= / (v(x)-v)W(x,v) 

(x>+ Vm ,...,x'+y m +...+ V2 )£X m - 1 x ' + 



X 





2 n " 2 (r - 


- 1 2/2 1 - ... - \y m \ 


II \ \n— 3 

- [x - x -y 2 - . . . -y m ) -v) 




n-1 L. 
• • • \ym\ 


n-1 


x — x 1 — y 2 - 


• • • • - y m - (r - \y 2 \ - ... - \y m \)v\ 


2n-4 



x [x(o,t-(x,v))(t -Si - \y 2 \ - ... - \y m \)E(x,x- (r - s x - \y 2 \ - ... - |y m |K 
x' + y m . . . + y 2 , . . . , x' + y m , x')k(x - (r - s x - \y 2 \ - . . . - \y rn \)v, v u v) 
xk(x' + y m + ... + y 2 ,v 2 ,vi) . . . k(x +y m + ...+ y i+1 ,v i+1 ,Vi) . . . 
k(x' + y m + y m -i, v m -i,v m - 2 )k(x' + y m , v', v m -i)S(x', v') 

\v(x') ■ v'\] x - x :- y ^..^y m ^ iT ^ s ,- ly , >l -...- ]ym[)v dy 2 . • .dy m dv, (8.13) 



s _ \x-x'-y2-----ym-(T-\y2\-...-\ym\)v\ 2 

1 2 ( t -\y2\-----\Vm-l\-( x - x ' ' -y2----Vm-l)- v ) 
v l = ym 

\ym\ 

Va = t~^t i i= 2...m — 1 

% ml 



for t G M. and a.e. (x, x') G dX x OX and for m > 3. 

We are ready to prove Theorem HSU We prove fl3~T7j) . (gHHD and f l3~T^|) . Let r G (0,T) 
and let x G <9X, x' G <9X and a; 7^ x'. Set t — |x — We first look for an upper bound on 
\j 2 (t, x, x')\. Using (I8.12p and the fact that a is a nonnegative function, we obtain 

|72(t,*,*')| < 2 n ~ 2 \\W\\ OD \\S\\ OD \\k\\ 2 O0 J 2 (r,x,x'), (8.14) 

where J 2 and S hn (r,x,x') are defined by (Q and (jOl . From (|8A4l) and fj8T3l- (j8"3|l it follows 
that there exists a real constant C such that 

|7 2 (r,x,x / )| < CWWISIWIAC sup J 2 (s,*,A whenn = 2, (8.15) 

(s,2,2')e(o,T}xi n xi n 

s> I z — z' I 

and 

rU ~ A ■al7a(r,a:,^)|<C||W|| 00 ||5|| 00 ||A ! ||2 0l when n = 3, (8.16) 



(^-k-^l)fi + in(^ 



and 

I /In— 2 

^|72(t,x,x')| <C||W|U|S|U|A;||L whenn>4. (8.17) 



T — \X — X\ 

Let m > 3. Using C I8. 131) we obtain 

\ lm (r,x,x')\ < 2 n - 2 \\W\\ O0 \\S\\ O0 \\k\\ZJ m (T,x,x'), (8.18) 

where 

t t <\ f ^(T(y),w(y),0) , , , e1 „v 
J m {T,x,x)= | r— I 1 T—[dy 2 ...dy rn _ ll (8.19) 

,(r,x,x') 1 2/2 1 • • • 12/m-ll 



m — 2,n I. 



and J2 (resp. £ m -2. )) is defined by (18.21) (resp. 118. lft ) and where y = (y 2 , . . . , y m -i), 

r(y) = r-\y 2 \-. . .-\y m ~i\, t {y) = \x-x'-y 2 -. . .-y m -i\ &ndw(y) = x-x'-y 2 - . . .-y m -i 
for y G (R n ) m - 2 . 
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Assume n = 2. Then using (I8.19p . (18. 3p and spherical coordinates (and (18.11) ). we obtain 
J m (T,x,x') < sup J 2 {s,z,z') / | — j — -: -dy 

( s ,z, z ')6(0,T) x i»xi" Je m - 2 n (r,x,x') 1 2/2 1 • • • |2/m-l| 



s> z — z' 



< (27r) m 2 sup J 2 (s,z,z') / ds 2 ■ ■ ■ ds m -i 

(s,z,z')G(0,T)xE™xR" J S 2 + ---+ S m-1< T 



s>|z-z'| 



s: >0, i — 2...m— 1 



= (27r) m - 2 / sup J 2 (s,2,0- (8-20) 

^771 ZJ! ( s ,z,z')e(0,T)xM" xi" 
s>|z— «'| 

Finally combining (18.201) and (18.181) . we obtain 

T m-2 

\lm(r,x,x')\ < (2ir)™-^W\U\S\U\k\\Z-, sup J 2 (s,z,z'), (8.21) 



^771 — ZJ! ( s ,z,z')6(i),T)x«"xl" 

s>|z-z'j 



Statement (13371) follows from (18TT5D . f[8~2l~l) . 

Assume ri > 3. Note that 



a — \z — z'\ \z — z'\ 

- 1 l - = 1 - J ^ < 1, 8.22 

and 

l / 1 / /il / 1 \ \ 2 / / 1 i \ \ 2 

I 2 ~ \ z ~ z I i 1 Z^+l- 2 - 2 I \ \ Wi / 1 + s 



1 + In ^ } < sup (1 - s) 1 + In , (8.23) 

V VM — k — ^1/7 se(o,i) V v 1 -^// 

for (/!, 2, z') G (0, T)xR"x E n such that |z - z'\ < \i. 

From fl87L9D . (ET221) and (E3D, fl8T23l) and (ESI), it follows that there exists a real constant 
C such that 

Jm(T, x, x') < CJ rn (T, x, x'), (8.24) 

where J m is defined by (18.61) . 

Assume n = 3. Combining (I8.18p . (18.241) . (18. 7p . we obtain that there exists a real constant 
C (which does not depend on r, x, x' and m) such that 



T — \x — X n 



hm(r,x,x')\ < CWWUSlUkWZ , ' _,, ' 1 + ln 



r + 


;r 


-x'\ 


r — 


,r 


— x'\ 



,n— 1 /A/Al fSu- 1 \ _\ m 3 



|x — a/| \ \r — |x — x'\) ) (m — 3)! 

(8.25) 

Statement (|3~T8j) follows from ffBTTHI) and (l8~25l) . 

Now assume n > 4. Combining (18.181) . (18.241) . (18. 8p . we obtain that there exists a real 
constant C (which does not depend on r, x, x' and m) such that 

| 7m (r, x, *')l < C"||W|U||^|U||A;||-(r - |x - s'|)|s - xf"™ 1 ( ^ 3) , ) ) . (8.26) 

Statement (I3TT91) follows from fl8~T7I) and (I8T26D . 

We now prove (I3.20p -( |3.2ip . Let n > 3 and m > 2. From the expression of 7 m (see (18. 12ft - 
dHU), it follows that 

h m (r,x,x')\ < 2 n - 2 \\W\U\S\U\k\\ZUr,x,x') (8.27) 



31 



where 



J mr ■■■ \y m \ n 



(y2<---<ym)ee m _ ln (T,x,x' ) 

(V +y m ,...,z'+E™ 2 «i)6(suPPfc) m_1 



where N and £ m ^x,n( T , x i x ') are defined by (17.21) and (18. ip . Note that 

|y m |>5 and r - |j/ f | > |x - x' - y 2 - . . . -y m | > 5, (8.29) 



i=2 

for (y 2 , ■■iVm) G ^m-i,n( r ) such that x' + y m G supp/c and x' + y 2 + ■ ■ ■ + y m G supp/c since 
suppA; C{zgI infygax |y — z\ > 5}. 

We prove fl£M . Assume n = 3. Using (1QH|) - (IQ9]1 . (TT4D and (l8TT0l we obtain 

J 2 (r, x, x) < 2vr5" 4 J B(r, x, x). (8.30) 

Therefore using (18. lip and ( 18. 30ft we obtain 

sup {s-\z-z\Y l (l + \n( S+ \*~ Z }X \ h(s,z,z') < oo. (8.31) 

(s,z,z')e(o,T)x8xxax \ \S 1 2 z \ J J 

s>|z-z'|>0 

Now assume m > 3. Using (18. 28ft . we obtain 



(y3,---,ym)eS m _ 2tn (T,x,x' ) 



J.(r,^< / J.(r-S^lwUx'+g f »,)d ft ....d«, | (g32) 



where J 2 is defined by ( 18.21) . Using (I8.32p and (18. 4p and the estimate sup rg( - 0jl ) r(l — ln(r)) 2 < oo 
we obtain 



I m (T,x,x') < D J 



(\x - x' - y 3 - . . . - y m \ 








" H r - T,7=3 \yi\)\ x - x ' -ys- ■■■-y m \ 



2 



(z'+E£L 3 J/i,z'+ym)6(suppfc) 2 

(8.33) 

where D := sup re(01) r(l-ln(r)) 2 sup (s , 2 , z /) G (o,T)xaxx9x (s— |z— z'|) _1 s|z— z'| (l + ln ( * ) 
J 2 (s, z'). If m = 3, then using (I8.29P with "(y 2 , . . . , y m )" replaced by "(1/3, . . . , y m ) n , we obtain 

J 3 (t, x, x') < 2t<T 4 j DVo1(£i i3 (t, x, x')) (8.34) 

(we also used the estimate |x — x' — y 3 \ + r — [ 3/3 1 < 2r for y 3 G £1,3(7", x, x')). If m > 4, then 
using (18.291) with "(y 2 , . . . , y m )" replaced by "(z/ 3 , . . . , y m )" , we obtain 



dy 3 ... dy. 



m—l 



I m (r,x,x') < 2r5- 4 DVo\(S h3 (r } x,x / )) [ 

y Is/ar • • • \ym-i\ n 1 

(S/3v.S/m-l)e£ m _3,3(T,a:^ / ) 

< 2rr 4 J DVol(^ 13 (r,x,x'))Vol(§ n - 1 ) m - 3 / ds 3 ...ds m _ 1 

t /(«3.-.«m-l)6(0,+oo)'»-3 



2rr 4 DVol(g"' 1 ) m " 3 ( J_ 3)! Vol(g 1 , 3 (r, x, x') (8.35) 
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3 + ... + s m _ 1 <T 

m— 3 



(we also used the estimate \x - x' - y 3 - . . . - y m \ +t - \y 3 \ - . . . - \y m \ < 2r for (y 3 , . . . , y m ) G 
£m-2,3( r ) %, %') and we performed the changes of variables yi = Sjd;,, (si,Ui) G (0, +00) x S n_1 ). 
Statement fl3~2UD follows from tKTTb . (IOT]) and ([SIMD (JHSSD (and <K^ ). 
We prove fl3T2TD . Let n > 4. Using (E2HD and (1731) . we obtain 



I m (T,X,x') < \\s\z-z'\ n 2 N(s,Z, z')\\L^(R s xdX z xR^) 

dy m ...dy 2 

\ ym \^\x -x'- YZ2 y*\ n ~ 2 (T - ET=2 \vi\) 



v / fym ■ ■ ■ ^y2 / s n(?\ 

I L, |n-l I-. I — V^m | T ,_9/_ v^m I.. IV \°-^) 



x' + y m esuppk 



Assume m = 2. Using (18.291) and (18.361) . we obtain 

I m (r,x,x') < 5- 2n+2 \\s\z ~ z'\ n - 2 N(s,z,z')\\ L o C(RsXdXzX ^^ (8.37) 
Therefore using (18. 9p . we obtain 



I m {r,x,x') < ^ 2 " +2 || S |z-zr^ 2 ^(5,^^IU-( R ,x e ^x R: ,)Vol n _ 2 (S n - 2 )7r(r+| a ;-x'|) f . 

(8.38) 

Assume m > 3. Using (18.291) and (18.36!) . we obtain 



/ TO (r,x,x') <6- 2n+2 \\s\z-z'\ n - 2 N(s,z 7 z')\\ Lac{Rs><dX2>< ^ i) 



dy m ■ ..dy 2 



[Xfal^-.-bm-lh- 1 " 

(j/2,---,J/m)e£'m-l,n(T,X,a;') 

(8.39) 

Note that |y m | + |x-x'-?/ m | < |y 2 | + - • . + \y m \ + \x-x' -y 2 - ■ ■ --y m \ for (y 2 , ■■■,y m ) G (M n ) m_1 . 
Hence 

ll/al +••• + bm-i| < t- bm| and |y m | + \x - x' - y m \ < t, (8.40) 

for (y 2 , y m ) G £ m _i, n (T, x, x') (see (EI])). Therefore 

lUr,x,x') <5- 2n+2 \\s\z- z'r 2 N(s,z,z')\\ L o, {RsXdXzXR n^ 

dy m ...dy 2 



y m es hn (r,x,x>) JTZt lwl<T-|wm| \V*\ n 1 • • • \Vm-i\ n 1 



„( W ') (m-2)! 



m-2 

<r M Vol(§"- 1 ) m -l S |z-^r~ 2 ^(s,z,/)|Uo 0(RsXaXzXM n ) Vol(f 1 ,4r,x,x / )) ' 



(m - 2)! 

r^VoKS"- 1 )™- 2 ^^ - z'\ n - 2 N(s, z, z')\\ Loo{RsXdXzXRni } Vol n _ 2 (§"- 2 )7r(r + |x - x'l) 



— - \ n_1 



(8.41) 



(m-2)!' 

Statement fl3~2TD follows from (15371) . flg^gj) and flHUj) . □ 
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9 Proof of Lemmas ED, EM lO and S3 



We remind the following change of variables for the proof of Lemmas 18.11 18.21 18.31 and I8.4L 

f(y)dy 



£l,n(T,t V,0)) 
( 



f 



t + s cos(<y?) J s 2 - tl . 
, — sin ip 



(0,27r)x(tQ,r) 



x il^£M dsd ^ ifn = 2> 



4a/s 2 - t 2 , 



f 



t + s cos((p) \/s 2 -t 2 



(9.1) 



sin (feu 



gn-2 x ( ,7r)x(t ,r) 



sin 



n-2 



t 2 COS 2 (</9) 



dudsdip, if n > 3, 



for / G L 1 (R n ) and (r, t , u) e (0, +oo) x (0, +oo) x S™" 1 such that r > t . 

Proof of Lemma\8J\ We prove (ET3jl . Let (//, z, z') G (0, T) x dX x ]R 2 be such that fi > 
From ([H2D and (J73D it follows that 



Z — Z . 



J 2 (n,z,z') 



2tt 



2tt 



dy 



i, a (M>(i,o),o) MvV" M) 2 - l*o(l,0) -y| 



(9.2) 



where t = | z — | . 

Using the chan 
s G (to, yu), we obtain 



Using the change of variables y = ^(1, 0) + (s cos(y?), v " 2 L ° sin (93)) (see (19.11) ). </? G (0, 27r), 



where 



2 

J 2 ,i(/i, s, ip)d(pdi 
s — t cos (93) 



J 2 ,l(/i,S,V?) 



(9.3) 
(9.4) 



a/s 2 - toV/ 1 - Sy/fM-to COs(^) ' 

for G (0, 27r) and s G (to,/i). 

We give an estimate on J 2 ,i- From (19.4D and the estimates /1 — tocos(y2) > s — t cos(ip) 
s + t > s — t cos(yj), it follows that 



J 2) l(/i,S,^) < 



(9.5) 



for <p G (0, |) and s G (to, A*)- Performing the change of variables s = t + e:(/i — t ) we have 



to 



1 



e(l-e) 



< +00, 



(9.6) 



for s G (t ,/i). Combining (19. 3p . (19. 5p . (19. 6p . we obtain 



sup 

((i,!,z')6(0,T)x3XxI 
fi~>\z—z f I 



J 2 (/i, 2;, z') < 87T Z 



^e(l - e) 



de < 00. 



(9.7) 
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Statement (JE2]) follows from (l9~7j) . 

We prove (IQjl . Let (/i, 2, z') G (0, T) x dX x ]R 3 be such that [i > \z - z'\. Set t = \z - z'\. 
From (E2D, (El and fl9TT]) . it follows that 

27rln f *H«l+l*-*'-«l 

/ I ^/.~"|iy|" — ■ I z z^ y I 



o^u, / At-|y|+|to(i,o,Q)-y| 
Z7rm U-|s/|-|io(i,o,o)-y| 



3 (M,to(l,0,0),0) \y\2(ti-\ y \)\t (i,o,a)-y\ 



dy 



£1 

2 



77 J2,i(l^,s } (f)dipds } (9.i 

Jt Jo 



where 

sin((^) In 

J 2> i(iJ l ,s,(p) = 



^— to cos(i/j) 



(s + to cos(ip))(2fi — s — 1 cos((p)) 



= Inf ^^^ V , Sin( ^ } ? +— ^ , A (9.9) 

\ /U — s / \2 / u(s + t c °s(y?)) 2/i(2/i — s — t cos(<£>)) / 

for <y9 G (0, 7r) and s G (to, /•*)• From (19.91) and the estimates 2ji — s — t cos(</>) > // — to cos((p), 

q < ln j Vfocos( y ) j < ln ^ it foUows that 



J 2)1 (/j,,s,ip) < In 



yU + t \ / sin(y9) sin(<£>) 



ji — s J \2/j,(s — t cos(y?)) 2/i(/i — t cos((/?)) 
for <p G (0, 7r) and s G (t ,fi). Therefore 

J 2l i(lJ,,s,ip)d(p < — ^ (\n ( S + t ° ^\ +\n f ^ + to 



2/it V V s — to J \fi — t 



'0 



ln 



M+*0 

fl—S 



t L + t o\ , n ( H + to 



< 1/ ^ ^ + in ^ , (9.10) 

We remind the following integral value 

" In ( E±±°\ ds = U- t ) In ( +fi-t . (9.11) 

to \V- S J \V-toJ 

Using the estimate ln f^±j) < ln (^r^ 1 ) for s G (t , ^), we obtain 

" In f ^ In ds = 2 A In ( ^] In f ^ rfs 



to \S-t J \n-sj J to \S-toJ \fJt- s 



^ ~ to J Jt \S — to 

^-<°>H^) +M2 0K^> + 1 >' (M2) 
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Combining (ra- (19~T2l and (1511]) . we obtain 

Statement follows from (15351) . 

We prove (1531) . Let n > 4. Let (/i, 2, z') G (0, T) x dX x R n be such that /i> \z-z'\. From 
(IO) and it follows that 

J 2 (aw') <<? / | 2/ | 1 -( / u-| z/ |)-^-/- z/ | 2 -^ = C / I?/! 1 — (^x— | 2 /|)- 1 |t (l, 0, 0)— ^l 2 —^, 

J E^ni^z') J£ 1 , n (^to(l,0...0),0) 

(9.14) 

where t = |* - z'\ and C = sup (A) j jF)e(0)T)xaXxffi „ /2|z - 5'| n - 2 iV(/i, 5, 5'). 

Using the change of variables y = y(l, 0) + (s cos(y), — 5 — - sin(yj)u;) (see (19.11) ). eg 6 (0, 7r), 
s G (to, A 4 ), ^ ^ § n ~ 2 , we obtain 

J 2 (a*, z, z') < C / / J2,i(n,s,<p)d(pds, (9.15) 
where C" = 2 n - 2 Vol n _ 2 (§ n - 2 )C and 



to ^0 



(s 2 -t 2 )^sin"- 2 (y9) 

2,l(// ' S ' ^ ~ (s + t cos(^))"- 2 (2/i - s - t cos(^))(s - t cos(^))"- 3 ' (9 ' 16) 

for ip G (0, 2tt) and s G (to,/i). 

We give estimates on J 2) i. Let y? G (0, |) and s G (to, A*)- From (19.161) and the estimates 

s + t cos(<£>) > s, a/s 2 — tg sin(<£>) < s — to cos(^) and the estimate 2fi — s — t cos(y?) > 
s — t cos(yj), it follows that 



s n 2 (s — t cos(y?)) 2 



s n x (s — t cosyep)) 

where Co is defined by ( 15.331) (we also used the estimate s — tocos(^) > s(l — cos( y))). Let 
<p G (| , 7r) and s G (to, A 4 ). From (19.161) and the estimates s — t cos(y?) > s, a/s 2 — t| sin(^) < 
s + to cos(y) and the estimate 2/i — s — t cos(^) > s, it follows that 



«M, ? )<J»<C,^P 71T , (9.18) 



s + to cos(<£>)) 2 s n : (s + t cos(<^)) ' 



where C is defined by ( 15.331) (we also used the estimate s + t cos(<p) > s(l + cos(y))). 
Combining (19~T7]) and (EOS]) and (171]) . we obtain 



Note that 



/ J 2 ,i{v,s,(p)dip < for s G (t , n). (9.19) 

Jo s 



— ^— c/s = — r n 2 ~ *° 9 ) = y /i - i -t +2 - 

< (9-20) 
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(we used the estimate t < \i which gives /i 1 1 £q +2 n < fi Hq 2 for i = ... n — 3). 
Combining (05]) . (IQfl - lEOOl) . we obtain 

J 2 (v,z,z')<C'^. (9.21) 

where C" does not depend on /i and z, z'. Statement (18.51) follows from (I9.2ip . □ 

Proof of Lemma\KM Let {r,x,x') E (0, T) xJETxlR™ such that r > \x-x'\ > 0. Set i = \x-x'\. 
Let (y 2 , • • -,y m -i) e ). Then t < |y 2 | + . . . + |j/ m _i| + |x - x' - y 2 - . • 

Therefore either |x — x' — y 2 — • • • — Um-il > z^ti or |x — x' — y 2 — . . . — y m -i | < -^-r an d there 
exists j G N, j = 1 . . . n such that \yj\ > z£br- Therefore using ( 18. 6ft we obtain 

1 " dy 2 ...dy m _i 



lit ± n 

J m (r,x,x') < ^2 /( 



|n-2 



i=2 - M-e-^K^ \y2\ n 1 ...\y m -l\ n X \X-x'-Y,i=2 Vi 
+ ^m,o( r ; x i x ) 

(m - 2) J m ,i(r, x, x') + J m , (r, x, x'), (9.22) 



where 



/ ' / dy 2 ...dy m _ 1 
J m0 (T, x,x) = / , i , (9.23) 



^i=2 m— 1 



f / /\ f dy 2 ... dy m _ x 

J I— '-E™2 »il<^T l^ 2 l •••|2/m-l| |X — X — 2^j= 2 

We first reduce the estimate of J m $ and J m> i to an estimate on 

Pm(r,x t x'):=j | ,S---f m -;_ r 0.25) 

From (19.23P and the estimate |x — x' — y 2 — . . . — y m -\\ > ^ry it follows that 

1 \ n-2 
TO — 1 



J mi o(r,x,x') < ( — — j P m (r,x,x'). (9.26) 

From (I9.24p and the estimates \y 2 \ > -^y > |x — x' — y 2 — . . . — y m -i| it follows that 

J i i\ s f d y 2 ' --dym-i , n ^ 

J mtl {T,X,X ) S /(»2.-,Vm-l)6e m -a,n('-.«.-') i — 1 1 — 7 ; ^rn-l uT~\ ' V*' Z '> 

l»2l>^T 

Therefore performing the change of variables "y 2 "= x — x' — y 2 — . . . — y m -i we obtain 

~ ~ ( 171 — 1 \ n 2 

J^iCr^x') <J m , < ( — J P m (r,x,x'). (9.28) 
Now we estimate Ps(r, x, x'). From (19.1 ft and ( I9.25P it follows that 
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Let n = 3. Then using the estimate cos (ip) > — 1 (and the fact that s , ^° cos ^ = s f ° + 

to(l-coB( y )) < j 2to J we obtain 

S + CQ COS(ip)) S+IQ COS(ip) ' 



P 3 (t,X,x') < 7T 



n(sm(<p) — 2 — ln(s + to cos(v?) j dipds 
V ' / 

27r I r — t + / hi(s + t )ds — / ln(s — t )ds J 
V Jt Jt J 

'r + t 



= 2n(r-t ) ^2 + \n^-—^jj (9.30) 

(we used the estimate ln(s + to) < l n ( r + ^o) f° r s £ (to? r ) an d we used the integral value 
( 19. lip ). Let n > 4. Using (19.291) and using the estimates a/s 2 — tg sin(<^) < s + tocos(y?), 
s + to cos(v?) > s(l + cos(</?)) and s — t cos(<y2) < s + to < 2s for (s, y?) G (t , r) x (0, tt), we 
obtain 

P 3 ( W <) < V^,(^) f r ^V)^( S -yos( P )) 

Jo J to 2(s + t cos(<^)) 2 

< CoVol n _ 2 (§"- 2 ) T r ~^ dyds < 7rC Vol ra _ 2 (S"- 2 )(T - t ), (9.31) 

A, A) s + t cos(v7) 

where Co is defined by (15.331) (we also used (17. ip ). 
Finally let m > 4 then 

(y 2 , • • • , Vm-i) e ^m-2,n(T, =^ I 2/2 G £i,„(r, x, x') and \] < r J . (9.32) 

Therefore using (I9.25j) . spherical coordinates ("j/i"= s^, Sj G (0, +oo), G § n_1 for i = 
3, . . . , m — 1) we obtain 



P m (r, x,x') < P 3 (r, x,x' 



%3 ■ • • 



|n-l 

_m— 3 



|M|+...+|]/m-l|<T 1 • • • bm-l| 

/_m—3 
n™" 1 ^ = P 3 (r,a;, a r / )Vol n _ 1 (S n - 1 ) wl - 3 7 — . (9.33) 
S3 +...+ Sm _i<T (m — 6)1 

0<Sj, i— 3, . . . ,m — 1 

Finally statement (JHUD follows from (l9~m ([9~2oD . (I9~25j) . (I9~3"0l and (ET351) . and statement 
(ED follows from (15351) . flQB) , (19~5B and (j935j) . □ 

Proof of Lemma \8. 31 Let n > 2. Using a rotation and (18. ip . we have 

Vol n (f ljn (r ) x, x')) = Vol n (£ lin (r, t„ei, 0)), (9.34) 

where t = |x — x'| and ei = (0, . . . , 0) G W 1 . 
From (I9.ip . it follows that 

vou^^toe.o)) = voi„_ 2( s"- 2 ) r r ( ^)f^ r a ' a -w;) ^. 



(9.35) 
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From (19.35)) and the estimate sin((p)^s 2 — 1\ < \J r 2 — t% for s G (to, T ), we obtain 



n-2 



vouvmo^o)) < voi„_ 2 (§«- 2 )(^-3) r r * 2 t ^^ dsd (p 

\ 2 / L J Q 4 x /s 2 - 1 2 



n-2 



< ivol n _ 2 (§ n - 2 ) ( y —^ | Vol(fi, 2 (r,toei,0)). (9.36) 



We remind that Vol (£^2(1", t ei, 0)) = 7r ^" +T - > v / ' r *o _ Therefore (18. 9p follows from (I9.36p . Lemma 



is proved. □ 
Proof 0/ Lemma \K4\ Let (//, 2, 2/) G (0, T) x OX x M 3 be such that /i > \z — z'\ > 0. Using the 



change of variables y = y(l,0) + (s cos(ip), — - sin(<^)u;) (see (19. II) ). <p G (0, 7r), s G (£0,^), 
w G S 1 , we obtain 

B(/i,z,z) = - / B 1 (n,s,tp)d(pds, (9.37) 

4 ./in ^0 



where 



^(/i, s, <p) = (s 2 - tl cos 2 (^)) sin(^) In ( ^ ) , ( 9 . 38 ) 



for (p G (0,2vr) and s G (t ,/i). Using f[9~3~g]) and the estimates In ( e=±jJ^M \ < i n (b±^l 
s 2 — tQCOs 2 (v?) < /i 2 , we obtain 

/ Bi(/2, s, (p)d(p < /i 2 / sin(y?)d(yg In | — - J, (9.39) 
Jo </ \ f 1 ~ s J 

for s G (t , M)- Combining f[9~3TD . flj£3gj) and (TDTTTT) we obtain 

/J. 2 tt / f fi + t 



B(ji, z, z') < ^(A* - *o) (h + 1 ) > ( 9 - 4 °) 

which proves (18.111) . □ 



10 The distributional kernel of the operators H m and the 



proof of Proposition 13.1 



Before we prove Proposition 13.11 we shall introduce and prove Proposition 110.11 given below, 
which gives the distributional kernel of the operators H m defined by ( 12.120 . 

Let E denotes the nonnegative mesurable function from W 1 x R" to R defined by 



E( Xl ,x 2 ) = e~ f ° X1 ^ a ( xl - a &=%l'&=& )ds G(xi,X2), for a.e. (x h x 2 ) 6l"x R n , (10.1) 

where B is defined by (I2.10p . For m > 3, we define recursively the nonnegative measurable real 
function E( ]_ i • • • i 0C 'i ll, 

) by the formula 

E(xx, ...,x m ) = E(x 1 , . . . ), (10.2) 

for (x 1 ,...,x m ) G (R n ) m . 

Concerning the distributional kernel of the H m , m > 2, we have the following result. 
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Proposition 10.1. We have 

H M M=I Pm(t, x, v, v')<p(x',v')dx'dv', (10.3) 

JxxS™- 1 

fort G (0,T) and a.e. (x,v) G X x S n_1 and for m > 2, where 

P 2 (t,x,v,x ,v) = / X(o,t~s 2 ){\x - {x - s 2 v )\) ■ 2n _ 4 



\x — S2V' — X' — {t — S 2 )v\ 

x \E(x, x — (t — Si — s 2 )v, x' + s 2 v' , x')k(x — (t — si — s 2 )v, Vi,v) 
xk(x' + s 2 v',v',vi)] _ x - S2V i- x '- {t - ai „ S2)v ds 2 , (10-4) 



n 



_ | X — 32 v f — x' — (t — 32 )v | 2 
^ 2(t — S2~ {x — x ! — S2V f )-v) 

fort G (0,T) and a.e. (x,v } x' } v') G X x § n_1 x X x § n_1 , and u>/iere 



P m (t,X,V,x',v') = / / X(0,t-s m -...- S2 )(k / + Sm^ + ... + S 2 ^2-^|) 

J/gn—lNm — 2 y s 2 + --- + s m<t 
* ' s^>0, i— 2...m 

2 n " 2 (t - s 2 - . . . - s m - (x - x' - s 2 v 2 - ... - s m _ii; m „i - s m v ') ■ v) n ~ 3 rs . , . 

X : ; ; ri^~4 [E(x,X-(t-S 1 -...-S m )v, 

\x-x' - s 2 v 2 - ... - S m __iW m _i - s m v' - (t - s 2 - . . . - s m )v \ 2n 4 

x' + s m v' + S m _iV m __i + . . . + s 2 v 2 , x' + s m v' + S m _iW m _i + . . . + S 3 V 3 , 

...,x' + s m v', x')k{x - (t - si - . . . - s m )v, vi,v)k(x' + s m v' + s m _it> m _i + . . . + s 2 v 2 , v 2 , v x ) 
■ ■ ■ k{x' + s m v' + s m _if m _i + . . . + s i+1 v i+ i,v i+ i, Vi) ... 

k{x' + s m v',v', f TO _i)] 3! -^- 32U2 -...- Sm . lt , m , 1 - Smt ,/-( t - 31 -...- 3m )„ ds 2 . . . ds m dv 2 . . . dv m -i, (10.5) 



_ |3:-3:'-S2U2----- s m-l"m-l- s m»' / -(t-S2-...-s m )u| 2 
1 "2(t-S2-...-s m -(x-x'-s 2 V2----s rn _ 1 v rn _ 1 -Sm.v')-v) 

forte (0,T) and a.e. (z, u, rr', «') G X x S™" 1 x X x S™" 1 , m > 3. 
Proof of Proposition \10.1\ Note that 

H 2 (t)(j)(x,v) = (J ^ £/i(t-Si-s 2 )A 2 i7i(si)A2i7i(s2)0tZs2dsiVar,«) 

C/i (t - Si - s 2 )A 2 £/i(si)A 2 cfsi J U 1 {s 2 )(j)ds 2 J (x, u) 
VJO / / 

i ft~S2 r- 

E(x, x — (t — si — s 2 )v) / k(x — (t — si — s 2 )v, f i, v) 

XE(x — (t — Si — S 2 )v, X — (t — Si — S 2 )v — SiVi) 

x / k(x — (t — s 2 — si)v — SiVi,v 2 , vi) 




^0 



xE(x — (t — Si — S 2 )v — SiVi, X — (t — Si — S 2 )v — S\V\ — S 2 V 2 ) 

X(f>(x — (t — si — s 2 )v — SiVi — s 2 v 2 , v 2 )dv 2 dvidsids 2 , 

for t G (0,T) and (x,v) G X x S"" 1 , where functions E are defined by ffT03j) - (TlT2]) . 
Using the change of variables u y(si,Vi) — (t — s 2 — Si)v + SiVi" we obtain 

H 2 (t)(j)(x,v) = [E(x, x - (t - si - s 2 )v,x - y,x - y - s 2 v 2 )k(x - (t - s x - s 2 )v,Vi,v) 

Jo Js™- 1 



x k(x - y,v 2 ,vi)} y- {t -s,- S 2)v 
1 s i 



_ \y-(t-s 2 )v | 2 
S1 ~ 2(t-s 2 -yv) 



x : ; : , 2n -4 <P{ X ~V~ s 2 v 2 ,v 2 )dydv 2 ds 2 . 

\y- {t-s 2 )v\ 
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Hence we obtain (110.31) . Note that 



(H 3 (t)<p)(x,v) = / H 2 (t-s 3 )A 2 U 1 (s 3 )<f ) ds 3 




Hence 



o 

t 

[3 2 {t - S 3 ,X,V,X 2 ,V 2 )(A 2 Ui(s 3 ))(/)(x 2 ,V 2 )dx 2 dv 2 ds 3 

JXxS"- 1 

(3 2 {t - s 3 ,x, v,x 2 ,v 2 ) / k(x 2 , v',V 2 )E(x 2 ,X 2 - s 3 v') 
XxS"- 1 Jo 7s™- 1 

x<p(x 2 — s 3 v ', v ')dv ' ds 3 dx2dv2- 



(H 3 (t)(j))(x,v) = / f3 3 {t,x,v,x\v')(j)(x',v')dx'dv\ (10.6) 



where 



f3 3 (t,X,V,x',v') = / / X(0,t~s 3 -s 2 )(W + s 3 v' - X + S 2 V 2 \) 

Jn™- 1 Jo Jo 

2™~ 2 (t — s 2 — s 3 - (x - s 2 v 2 —x 1 - s 3 v') ■ v) n ~ 3 
\x — x' — s 2 v 2 — s 3 v' — (t — s 2 — s 3 )v\ 2n ~ 4 
x [E(x, x — (t — si — s 2 — s 3 )v, x' + s 2 v 2 + s 3 v', x' + s 3 v', x') 
x k(x - (t — s 3 — s 2 - si)v, Di, v)k(x' + s 2 v 2 + s 3 v' , V 2 , V i) 
x k(x' + s 3 v', v', v 2 )] ^ x - x i- S2 v 2 -s a vi-(t- sl -s 2 -s a )v ds 2 ds 3 dv 2 . (10.7) 

_ \x-x' -s 2 v 2 -s 3 v' ' -(t-S2-s a )v\ 2 
2(t-s 2 -s 3 -(x-x'-s 2 v 2 -s 3 v')-v) 

The proof of (110.51) follows by induction from (110.61) and (I2.13p . □ 
Proof of (18321) - (Km . We recall that 

{A 2 G-{s)(ps) (z,w) = ( [k(z,v',w)S{x',v')\u(x') ■ v'W vl _ z- x > , E ^ Z, ?J (/){s-\z-x'\, x')dfi(x'), 

JdX \Z — X \ 

(10.8) 

for a.e. (z, w) G X x S n 1 and G -^((0, rf) x <9X) (see the derivation of (13. 9p and f)3.10p given 
in Section [3]). 

Let m = 2. Then from fl2~T2l and (EJJ) it follows that 
A 2 ,s,w(<P)(t,x) = / (u(x) ■ v) W(x, v) / / / / [k(x - (t - s - sx)v,V!,v) 

+ 1 J-oo Jo JS"- 1 JdX 

X k(x — (t — S — Si)v - SiVi, V , V!)S (x , v')\is(x') ■ v'\] v , = x-^-s-s^v-s^-x' E(x,x — (t — s — si)v, 

\x — (t — s — s-±)v — s^vi — x f \ 

,,0(S- \X- (t - S - Sx)v - SxVx -X'\,X') AAA 

x — [t — s — Sx)v — S\V\, x ) j -. r 1 — : du(x )dvidsidsdv. 

\X — (t — S — Si)V — SiVi — x'\ n ~ L 

Performing the change of variables y(s\, vi) = (t — s — si)v + s\Vi, we obtain 
A 2 , s ,w((p)(t,x) = {v{x) ■ v) W(x, v) / X(o,t-s)(\y\) 

Js^xdXxM™ J~oo 

x [E(x, x — (t — s — s\)v, x — y, x')k(x — (t — s — si)v, vi, v)k(x — y, v', v i)S(x', v') 

I / /\ /,, 2 n ~ 2 (t -s-y v) n - 3 (f)(s -\x-y-x'lx') , inm 

x Hx ).v\] n= K-^ \{t-s)v-y\^\x-y-x>\^ ^ (1 °' 9) 

Vi _ y-( t -j>- s -i) v 

i x — ii — x r 
v =1 77 
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Performing the change of variables "y" = x — x' — y and t' = s — \y\ we obtain (I8.12p . 
Let m = 3. Then from (1331) . (ITTHj]) and (fPDTgl) it follows that 

A 3AW ((f>)(t,x)= f (v(x)-v)W(x,v) [ [ fo(t-3,x,v,x 2 ,v 2 ) (10.10) 

JS™" 1 J-oo JXxS"" 1 

E(x 2 x') 

[k(x 2 , v' ', v 2 ) S(x' ', u')|i/(a/) ■ t/|] ,_ xg-g' i ti -; 0(g - 1^2 - x'\,x')dfi(x')dx 2 dv 2 dsdv, 

ax v ~ \ x 2~ x '\ F2 - a? T 

for t e (0, T) and ac G <9X. From ffTIUOl and (117141) we obtain 

A 3yS ,w{<P){t,x) = (u(x) ■ V)W(X,V) / / / X(0,t-s-s 2 )(|#2 - (z - s 2 u 2 )|) 



n-1 



XxS"- 1 xdX J-oo JO 
n— 3 



2 n-2 ^ _ s _ g 2 _ ( x _ S2 ^ 2 _ . v } 

\x 2 — X'| n_1 |x — S 2 V 2 — X 2 — (t — S — S 2 )t> | 2 ™~ 4 

x x - (t - si - s 2 )v, x 2 + s 2 t> 2 , x 2 , x')k(x — (t — si — s 2 )v, Vi, v)k(x 2 + s 2 v 2 , v 2 , Vi) 
k(x 2 , v', v 2 )S(x', v')\v{x') ■ v'\) = x -s 2 v 2 - X2 -(t-s-s 1 -s 2 )v <f)(s — \x 2 — x'\,x')ds 2 dsdx 2 dv 2 dfi(x')dv. 

V l- Bl 

s _ \x-s 2 v 2 -x 2 -{t-s-s 2 )v\' 2 

1 2 ( t - s - s 2-( x - x 2- a 2 v 2)- v ) 
v , = x 2 -x> 
\x 2 -x'\ 

Performing the change of variables y 2 = s 2 v 2 and y$ = x 2 — x' we obtain (18.131) for "m = 3" . 
Let m > 3. From (TT0~3l . (TT031) and flOl it follows that 

A m +i,s,w(<j))(t } x) = / (v(x) ■ v) W(x, v) / / /? m (t - £' - |x m - x'|,x,t;,x m ,f m ) 

J$x~+ JdX J X J (— oo,t— |x m — x'|)x§ n_1 

[fc(x m , v',v m )S(x', v')\v(x ) ■ v'W r ,_ x m -x> i E ^ Xm ^J 4>(t' , x')dn(x')dt'dx m dv m dv 

7m+i(^ _ x')<p(t ', x')dt'dfi(x'), (10.11) 

(0,r;)xax 

where 

r/ m +i(T,x,x') := J ^ (u(x) ■ v)W(x,v) [ X(o,+oo)(r-\x m -x'\) 

x,+ 

/<? n —l\m — 2 J s 2 + ... + s m <T — \xm- 
y > ' Sj>0, i=2. ..m 



'XxS n - 

, X(0,r-\x m -x'\- Sm -...-s 2 )(\x m + S m V m + ... + S 2 V 2 - x\) 



2 n - 2 ( 


r — 


| 




- s 2 v 2 - 


- . . . - s m v m ) ■ v) n 


-3 




n-1 


% %m ~ 


- s 2 t>2 - ... - s m w m - (r - | 




- s 2 - ... - s m )v\ 


2n-4 



x [E(x, x - (r - \x m - x'\ - Si - . . . - s m )v, x m + s m v m ... + s 2 v 2 , x m + s m v m + . . . + s 3 v 3 , 
x m + s m v m , x m , x')k(x - (t - \x m - x'\ - si - . . . - s m )v, v ± , v)k{x m + s m v m + ... + s 2 v 2 , v 2 , v x ) 
. . . k(x m + s m v m + . . . + s i+ iv i+ i,v i+1 , Vi) . . . 



k(x m + s m v m7 v m ,v m -i)k(x m} v' \v m )S(x' \v')\v(x') ■ v'\ 



_ X — X m -S 2 V 2 -... — 3 m V rn -(T—\x rn — x'\—3 1 —... — S m )v 
"1- Bl 

_ \ x-x m — s 2 v 2 — ... — s m v m -(T—\x m -x'\—s 2 — ...-s m )v\ 2 
1 2(t-s 2 -...-s m -(i-i'-S2»2--.Sm»m)'») 
L ./_ x m -x' 
\x m — x'\ 



ds 2 . . . ds m dv 2 . . . dv m -idx m dv m dv . (10.12) 

Performing the change of variables yi = Sjfj, i = 2 . . .m, and y m +i — Xm+i ~ x', we obtain 
(E33D for "m > 4" . □ 
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